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Abstract. The progress of the last decade in perturbative quantum field theory 
at high temperature and density made possible by the use of effective field theories 
and hard-thermal/dense-loop resummations in ultrarelativistic gauge theories 
is reviewed. The relevant methods are discussed in field theoretical models 
from simple scalar theories to non-Abelian gauge theories including gravity. In 
the simpler models, the aim is to give a pedagogical account of some of the 
relevant problems and their resolution, while in the more complicated but also 
more interesting models such as quantum chromodynamics, a summary of the 
results obtained so far are given together with references to a few most recent 
developments and open problems. 
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1. Introduction 

Ultrarelativistically hot and dense matter needs to be understood in many problems 
ranging from early-universe cosmology [1, 2] as well as the astrophysics of compact 
stars [3] to the current and prospected experimental programs at relativistic heavy- 
ion colliders in the USA (RHIC at Brookhaven) and in Europe (LHC at CERN), 
which seek to produce and investigate the new state of matter called the quark-gluon 
plasma [4] . The theoretical framework for analysing equilibrium and near-equilibrium 
properties is thermal field theory, which at a fundamental level involves Abelian and 
non-Abelian gauge theories. 

In this review we shall focus on the technical difficulties to derive results by 
analytical, weak-coupling techniques. Because of asymptotic freedom of non-Abelian 
gauge theories, this is a valid approach at least at sufficiently high temperatures and 
densities. Much of the present activities is concerned with quantum chromodynamics 
(QCD), which is however strongly interacting at all temperatures of experimental 
interest. The perturbative calculations that have been carried out to high loop 
orders for thermodynamic potentials seem to signal that perturbative QCD at finite 
temperature is loosing any predictive power below ridiculously high temperatures 
<~ 10 5 GcV, where the coupling is so small that everything is described well enough 
by an ideal gas of quarks and gluons without much need of perturbative refinements. 
However, recent studies indicate that this conclusion is much too pessimistic, and it 
appears that even in strongly coupled QCD it is indeed possible (at least to some 
extent) to employ weak-coupling methods if proper use is made of effective field 
theories and the fact that the nontrivial spectral data of quasiparticles can absorb 
much of the strong interactions [5]. 

In the case of static quantities at high temperature and small chemical 
potential the relevant techniques are effective field theories which take advantage 
of the phenomenon of dimensional reduction; for dynamical quantities the required 
techniques involve the effective theory produced by hard thermal loops [6]. At 
high chemical potential and low temperature, the resummation of hard dense loops 
becomes important even in static quantities and also for describing the astrophysically 
interesting case of colour superconductivity [7]. 

The present review may to some extent be viewed as an extension of the material 
presented in the textbooks of Kapusta [8], Das [9] and LeBellac [10], with a certain 
amount of overlap with the latter which already covers the most important aspects of 
hard thermal loop rcsummations. Of course, numerous interesting new developments 
could only be touched upon, but hopefully enough references are provided as pointers 
to the recent literature; a lot of topics had to be omitted by lack of space or competence. 

1.1. Outline 

In section 2, after a brief recapitulation of the fundamental formulae of quantum 
statistical mechanics, the imaginary-time and real-time formalisms in thermal field 
theory are discussed with an overview of some more recently developed alternative 
approaches in the real-time formalism and in the treatment of gauge theories. 

In section 3 we then exhibit some of the particular issues that arise in finite- 
temperature field theory in a simple scalar model: the appearance of thermal 
masses, the phenomenon of dimensional reduction at high temperatures, the need 
for a resummation of naive perturbation theory, and the problem of a deterioration 
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of apparent convergence after resummation and possibilities for its improvement. 
After a short discussion of the phenomenon of restoration of spontaneously broken 
symmetries at high temperatures at the end of section 3, section 4 describes the 
progress made during the last decade in the (resummed) perturbative evaluation 
of the thermodynamic potential of unbroken nonabelian gauge theory. At high 
temperature and not too high chemical potential, dimensional reduction provides 
the most efficient means for calculating pcrturbatively the thermodynamic potential, 
and this calculation has been carried to its limits, which are given by the inherent 
nonperturbative nature of nonabelian magnetostatic fields, ft is argued that initial 
pessimism regarding the utility of weak-coupling methods in hot QCD is unnecessary, 
and that the results rather indicate that perturbative methods, when properly 
improved, may work already at temperatures that are only a few times higher than 
the deconfinement transition temperature. The recent progress made for the case 
of a nonvanishing fermion chemical potential is also reviewed, which includes quark 
number susceptibilities and non-Fermi-liquid contributions to the low-temperature 
entropy and specific heat of QCD and QED. 

Section 5 discusses the structure of the propagators of a gauge theory with 
fermions, the leading-order results for the respective quasiparticles and the issue of 
gauge (in)dependence, both in unbroken nonabelian gauge theory and in the presence 
of colour superconductivity. 

Section 6 then introduces the concept of hard-thermal/dense-loop resummation 
which is typically necessary (though not always sufficient) to calculate the effects 
of collective phenomena such as dynamical screening and propagating plasmons. 
Two recent approaches to improve the problem of the poor apparent convergence 
observed in the thermodynamic potential are discussed, which suggests that already 
the lowest order calculations in terms of quasiparticles may capture the most important 
contributions even in the strongly coupled quark-gluon plasma at temperatures that 
are only a few times higher than the deconfinement temperature. 

The known results concerning next-to-leading order corrections to the quasi- 
particle spectrum in gauge theories are then reviewed in section 7. In this section 
a few simpler cases such as the leading infrared-sensitive contributions to the Debye 
mass and to damping rates and dynamical screening lengths are shown in more detail. 
Other cases of special interest are enhancements from collinear physics, which include 
non-Fcrmi-liquid corrections to the fermion self-energy in the vicinity of the Fermi 
surface and a modification of longitudinal plasmons for nearly light-like momenta. 

Section 8 briefly discusses some of the progress made recently in the case of 
ultrasoft scales, where nonabelian gauge fields have nonperturbative dynamics, and 
in the case of observables that are sensitive to collinear physics, where also infinitely 
many loops contribute even after hard-thermal-loop resummation. 

Section 9 finally considers thermal field theory in a curved background geometry 
and in particular the hard thermal loop contribution to the gravitational polarization 
tensor. The latter encodes the physics of cosmological perturbations in the presence 
of nearly collisionless ultrarelativistic matter, for which analytical solutions can be 
obtained in the physically relevant case of a conformally flat geometry. 

2. Basics 

In a relativistic quantum theory where interactions typically imply the destruction or 
creation of particles it is appropriate to formulate a statistical description by means 
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of the grand canonical ensemble. 

A system in thcrmodynamical equilibrium for which only mean values of energy 
and any conserved charges are prescribed is characterized by a density matrix p such 
that [p, H ] = with H the Hamiltonian operator and the requirement of maximal 
entropy 

S = ( — In p) = — Tr p In p. (2.1) 

The temperature T and the chemical potentials pi appear as Lagrange multipliers 
j3 = T -1 and oti — —f3pi determining the mean energy (H) and mean charges (Ni), 
respectively, with [H, Ni] = = [N i} Nj]. This singles out 

i 5 = Z- 1 cxp{-/3 J ff-^a 4 iVj, (2.2) 

i 

where the normalization factor Z is the grand canonical partition function 

Z(V,/3,pi) = Trexp{-/?iJ -^OjJVi}. (2.3) 

i 

The partition function, Z, determines all of the other conventional thermody- 
namic (or rather thermo-static) quantities such as pressure, entropy, energy, and 
charge densities, which are denoted by P, S, £, and A/", respectively. In the ther- 
modynamic infinite-volume limit (V — > oo) these are given by 

dP 

S ^ S/V= dT> (2 ' 5) 
B U V 1^. (2.6) 

Combining (2.1) and (2.2) one obtains the Gibbs-Duhem relation in the form 
S = (- hip) = 1 In Z + (3£ + ^iXi = P(P + £ - V^i) ( 2 ' 8 ) 

i i 

or 

E = -PV + TS + fM N i , (2.9) 

which explains why P was introduced as the (thermodynamic) pressure. A priori, 
the hydrodynamic pressure, which is defined through the spatial components of the 
energy-momentum tensor through |(T M ), is a separate object. In equilibrium, it 
can be identified with the thermodynamic one through scaling arguments [11], which 
however do not allow for the possibility of scale (or "trace") anomalies that occur in 
all quantum field theories with non-zero /3-function (such as QCD). In [12] it has been 
shown recently that the very presence of the trace anomaly can be used to prove the 
equivalence of the two in equilibrium. 

All these formulae pertain to the rest frame of the heat bath. A manifestly 
covariant formulation can be obtained by explicitly introducing the 4- velocity vector 
tt M of this rest frame. In thermal equilibrium the energy momentum tensor takes the 
form 

T" I ' = &V + P(uV-))'"'), (2.10) 



Thermal field theory 



6 



where rjV = diag(+, — , — , — ) is the Minkowski metric. Introducing the 4-vectors 
(i^ = (3u^ and j? = riiU^, the grand canonical partition function can be written 
covariantly as [13] 

Z = Trcxp f dS M {- J 9„f - a di} ( 2 - n ) 
J i 

with reference to a hypersurface X with normal vector u. 

It is of course more convenient to stick to the rest frame of the heat bath. Lorcntz 
invariance appears to be broken then, but it is usually a trivial matter to switch to 
Lorentz covariant expressions using u M . 

Information about the dynamics of a thermal system can be obtained by 
considering thermal expectation values of the generally time-dependent observables 
(in the Heisenberg picture). In linear response theory, the time evolution of small 
disturbances of an equilibrium system is determined by the correlation functions of 
pairs of observables [14]. 

For both the purpose of calculating the partition function Z and correlation 
functions there exist two equivalent but rather differently looking formalisms to set up 
perturbation theory. These correspond to the two most popular choices of a complex 
time path in the path integral formula 

(T c <pi • • • <p n ) = J\f J VfLpi ■ ■ - fn expi J dt j d 3 x£ , (2-12) 

c 

where T c denotes contour ordering along the complex time path C from to to to — i/3 
such that U e C, and t\ y ti > ■ ■ ■ > t n with respect to a monotonically increasing 
contour parameter. C is the Lagrangian, which in the presence of a chemical potential 
jU 7^ may be replaced by C — > C = C + fiAf, provided Af does not contain time 
derivatives. Analyticity requires that along the complex time path the imaginary part 
of t is monotonically decreasing [11]; in the limiting case of a constant imaginary part 
oit along (parts of) the contour, distributional quantities (generalized functions) arise. 

Perturbation theory is set up in the usual fashion. Using the interaction-picture 
representation one can derive 

(T c Oi ■ • • O n ) = ^(T c Oi ...O n e'fc £ ') , (2.13) 

where £7 is the interaction part of £, and the correlators on the right-hand-side can 
be evaluated by a Wick(-Bloch-DeDominicis) theorem: 

{Tce iJ c * xjv)o =e K ? {-\jJ c d i xd*x'j{x)D c {x-x l )j{x')}, (2.14) 

where D c is the 2-point function and this is the only building block of Feynman graphs 
with an explicit T and fj, dependence. It satisfies the KMS (Kubo-Martin-Schwinger) 
condition 

D c (t-if3) = ±c-^D c (t), (2.15) 

stating that e lflt D(t) is periodic (anti-periodic) for bosons (fermions). 

A KMS condition can be formulated for all correlation functions in thermal 
equilibrium, and can in turn be viewed as a general criterion for equilibrium [15]. 
There exists also a relativistic version of the KMS condition [16, 17], which encodes 
the stronger analyticity requirements of relativistic quantum field theories. The 
relativistic spectrum condition H > |P|, where P is the total three- momentum, 
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implies analyticity involving all space-time variables x —> z e C 4 in tube domains 
|Imz| < Imz < (3 — |Imz|, whereas the usual KMS condition corresponds to 
restricting this to Imz = 0. 



2.1. Imaginary-time (Matsubara) formalism 

The simplest possibility for choosing the complex time path in (2.12) is the straight 
line from t to t — i/3, which is named after Matsubara [18] who first formulated 
perturbation theory based on this contour. It is also referred to as imaginary-time 
formalism (ITF), because for to = one is exclusively dealing with imaginary times. 

Because of the (quasi-)periodicity (2.15), the propagator is given by a Fourier 
series 

D c (t) = — -Y^D{z v )e- iz »\ D{z u )= / dtD c (t)e iz ^ (2.16) 
iP v Jo 

with discrete complex (Matsubara) frequencies 

n ■ m f Z bos. ,„ 

z v = 2ttw/(3 + n, v e | z _ i {cim (2.17) 

Since D{z v ) is defined only for a discrete set of complex number, the analytic 
continuation to arbitrary frequencies is unique only when one requires that |-D(z)| — > 
for \z\ — > oo and that D{z) is analytic off the real axis [19]. Then the analytic 
continuation is provided by the spectral representation 

K ' 2ir k - z V ' 

The transition to Fourier space turns the integrands of Feynman diagrams from 
convolutions to products as usually, with the difference that there is no longer 
an integral but a discrete sum over the frequencies, and compared to standard 
momentum-space Feynman rules one has 

/ i^I^^E / (03< W*{k) - mirfS,o5 3 (*). (2.19) 

However, all Green functions that one can calculate in this formalism are initially 
defined only for times on C, so all time arguments have the same real part. The 
analytic continuation to several different times on the real axis is, however, frequently 
a highly involved task [11], so that it can be advantageous to use a formalism that 
supports real time arguments from the start. 



2.2. Real-time (Keldysh) formalism 

In the so-called real-time formalism(s), the complex time path C is chosen such as 
to include the real-time axis from an initial time to to a final time tf. This requires 
further parts of the contour to run backward in real time [20, 21] and to end up at 
to — i/3- There are a couple of paths C that have been proposed in the literature. The 
oldest one due to Keldysh [22] is shown in figure 1 , where the first part of the contour 
Ci is on the real axis, and a second part runs from tf — iS to t — iS with 5^0. 
For some time, the more symmetric choice where the backward-running contour is 
placed such as to have imaginary part —S = —(3/2 with two vertical contour pieces 
at t an d t/ of equal length has enjoyed some popularity [23, 24, 11] in particular in 
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Figure 1. Complex time path in the Schwingcr-Keldysh real-time formalism 

the axiomatic thermo-field-dynamics (TFD) operator formalism [25, 26, 27, 28] but by 
now the original Keldysh contour seems to be the one most widely employed, above all 
because of its close relationship to the so-called closed time-path formalism [22, 29, 30] 
of noncquilibrium thermodynamics. 

With the Keldysh contour, if none of the field operators in (2.12) has time 
argument on C\ or C2, the contributions from these parts of the contour simply cancel 
and one is back to the ITF. On the other hand, if all operators have finite real time 
arguments, the contribution from contour C3 decouples. The standard argument for 
this assertion [23, 24, 11] relies on the limit to — > —00 and the fact that the propagator 
connecting contour C\ and C 3 decays by virtue of the Riemann-Lebesgue theorem 



duj e - iw (*-*°) -^—— p (uj, k) ( 2.20 ) 



for A e (0, (i). However, as argued by Niegawa[31] there are cases where this line 
of reasoning breaks down. The decoupling of the vertical part(s) of the contour in 
RTF does however take place provided the statistical distribution function in the 
free RTF propagator defined below in (2.23) does have as its argument |fco| and not 
the seemingly equivalent ojk = (k 2 + to 2 ) 1 / 2 . Reference [31] also relied on the limit 
to — > —00, but Gclis [32, 33] pointed out that the decoupling of C3 should take place 
regardless of the magnitude of to, and showed this indeed to be the case provided the 
"n(|fco|) prescription" is used (see also [34, 35]). 

With only C\ and C2 contributing, the action in the path integral decomposes 
according to 

//'OO /*oo 
C(ip)= dt£{<pW)- dtC{<pW) (2.21) 
J— 00 J— 00 

C1UC2 

where one has to distinguish between fields of type 1 (those from contour Ci) and 
of type 2 (those from contour Ci) because of the prescription of contour ordering in 
(2.12)4 From (2.21) it follows that type-1 fields have vertices only among themselves, 
and the same holds true for the typc-2 fields. However, the two types of fields are 
coupled through the propagator, which is a 2x 2 matrix with non- vanishing off-diagonal 
elements: 

, _ ( (T^(tMt')) a(ip(t')<p(t)) \ 

( ^ M t')) (iv(tMt'))J- (2 - 22) 

Here T denotes anti-time-ordering for the 2-2 propagator and a is a sign which is 
positive for bosons and negative for fermions. The off-diagonal elements do not need 

% Type-2 fields are sometimes called "thermal ghosts", which mislcadingly suggests that type-1 fields 
are physical and type-2 fields unphysical. In fact, they differ only with respect to the time-ordering 
prescriptions they give rise to. 
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a time-ordering symbol because type-2 is by definition always later (on the contour) 
than type-1. 

In particular, for a massive scalar field one obtains 

»•<*>- (*£fe> M ^" 2) ) + ^ 2 -»><i^i)(! ;)• 

= D[,(fc) +D^(fc) (2.23) 

where n{w) = [e^ - l]" 1 and S ± (k 2 - to 2 ) = 6(±k )S(k 2 - to 2 ). The specifically 
thermal contribution is a homogeneous Green function, as it is proportional 
to S(k 2 — to 2 ). Physically, this part corresponds to Bose- Einstein-distributed, real 
particles on mass-shell. 

The matrix propagator (2.23) can also be written in a diagonalized form 
[36, 28, 37] 

V c (k)=M(k a )( [G n F ° Wo) (2.24) 



-iff t 

with G F = l/(fc 2 - m 2 + ie) and 

If e kP\ko\ e (S-\p)k \ 

MW = VTO(e(V* ) ^ 

where 8 — > for the Keldysh contour. 

For complex time paths with 8 <x (3, the T — > limit decouples type-1 and type-2 
fields completely. For <5 — > 0, the limit T — > leads to 

M(M^°M (fe )=(^ o) 6{ - ko) ) (2.26) 

so that one still has propagators connecting fields of different type. However, if all 
the external lines of a diagram are of the same type, then also all the internal lines 
are, because J\i 0(k9 i -.) = when ^\ k^ = and any connected region of the other 
field-type leads to a factor of zero. 

The matrix structure (2.24) also applies to the full propagator, and consequently 
the self energy ill = D 1 — Dq" 1 has the analogous form with M _1 in place of M, 

n(k) = M(ko)( U Q F _°^M(fc ). (2.27) 

It is also possible to diagonalizc in terms of retarded and advanced quantities 
according to 

D c (fc)=U(fc )( i( ;* i( 2.) V fo) ( 2 - 28 ) 

with 

^•>=0 -a:'X)))< v<M=( i+ f°> ••*>) (2.29) 

and to include the matrices U and V in the vertices [38, 39]. This has the advantage 
of leading to n-point Green functions with well-defined causal properties, which 
correspond directly to the various analytic continuations of ITF Green functions. On 
the other hand, the type-l/type-2 basis leads to rather involved relations [40, 41, 42]. 
Because of 

U(fc)n = V T (-k), n= (J jY (2.30) 



Thermal field theory 



10 



more symmetric retarded/advanced Feynman rules can be formulated by including a 
factor T\ in U and putting Gr and Ga in the off-diagonal entries of (2.28) [43, 44]. 

A precursor of this transformation is in fact given by the so-called Keldysh basis 
[22, 29, 44] 

<P+ = §(<Pi + <P- = (<Pi - V2) (2.31) 

(sometimes labelled by indices r, a instead) . This also has the advantage of a rather 
direct relationship to retarded/advanced n-point Green functions, and because the 
transformation does not involve n(fco), it is of use also in the nonequilibrium closed- 
time- path formalism [45]. 

Another economical method to derive retarded/advanced quantities in the real 
time formalism is provided by the use of the outer products of 2-component column 
vectors [46, 47, 48] as worked out in [49, 50]. 

2.3. Extension to gauge theories 

In the partition function (2.3) and in thermal averages (Q) = TrpQ, the trace is 
taken over the physical Hilbert space. But covariant formulations of gauge theories 
are built in larger spaces containing unphysical states, while the definitions following 
(2.3) are true only in the physical Hilbert subspace. The standard solution is to 
extend the trace to the larger unphysical space and to cancel unphysical contributions 
by Faddeev-Popov ghosts. 

In the path integral formalism, the Faddeev-Popov ghost fields arise from a 
functional determinant in the configuration space of the bosonic gauge fields. This 
requires that although Faddeev-Popov fields behave as fermions with respect to the 
diagrammatical combinatorics, they are subject to the same boundary conditions as 
the gauge bosons and therefore have the same statistical distribution functions, namely 
Bose-Einstein ones [51]. 

In operator language which starts from a BRS invariant theory involving fcrmionic 
Faddeev-Popov field operators [52], this prescription can be understood through the 
observation [53] that the operator expiTrN c , where N c is the ghost-number operator, 
is equivalent, under the trace, to a projection operator onto the physical Hilbert space. 
This means that the fermionic Faddeev-Popov fields are given an imaginary chemical 
potential \i c = in//3. But a Fermi-Dirac distribution with such a chemical potential is 
nothing other than a Bose-Einstein distribution. 

There is however an alternative approach, developed in [54, 55], which avoids 
assigning thermal distributions to unphysical degrees of freedom altogether. In the 
real-time formalism, one may switch off the gauge coupling adiabatically as the 
beginning of the time contour is moved to — oo. Then the condition to choose physical 
states can be the same as in Abelian gauge theory. The unphysical states are identified 
as those which are due to the Faddeev-Popov fields and the temporal and longitudinal 
polarizations of the gauge fields. Because the free Hamiltonian is a sum of commuting 
parts containing respectively only physical and unphysical operators, and because the 
unphysical part has zero eigenvalue on the physical states, all unphysical contributions 
factor out such that only the transverse polarizations of the gauge fields acquire 
thermal parts in their propagators. 

In Feynman gauge for instance the gauge propagator, which usually is simply 
T>^ u — —n^T) , with D c the matrix propagator (2.23), becomes 

D"" = -A" V D C - « I/ - A" v )T> c a (2.32) 
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with 

= 0, A ij = -(5 ij - — £-) (2.33) 

and the Faddeev-Popov ghost propagator remains non-thermal, Dpp = D§. 

In a general linear gauge with a quadratic gauge breaking term in momentum 
space according to 

4-br. = ~A^{-k)Uf v A v {k) (2.34) 
the vacuum piece generalizes by replacing 

rT ~» if" - — +Cf + ak 2 )^—. (2.35) 

f-k U \f-kf 

The ghost propagator is replaced by (/ • fc) _1 with real-time propagator matrix 
analogous to the vacuum part of the gauge propagator. 

At finite temperature, where manifest Lorentz invariance has been lost anyway, 
the modification (2.32) does not introduce additional non-covariance. In fact, it 
simplifies calculations of thermal contributions in general gauges [54] , but it makes it 
more intricate to investigate resummation effects [55]. 

In gauges where the ghost degrees of freedom are non-thermal anyway, such as 
Coulomb gauge or axial gauges, the above Feynman rules are identical to those of the 
conventional approach. In particular, it reproduces the real-time Feynman rules for 
temporal axial gauge of [56] which presents major difficulties in the imaginary-time 
formalism [57]. 

3. Resummation issues in scalar </> 4 -theory 

Before discussing gauge theories further, we shall consider perturbation theory at finite 
temperature in a scalar field theory with quartic coupling and address the necessity 
for resummations of the perturbative series. 

3.1. Daisy and foam resummation 

A particularly simple "solvable" model is given by the the large N limit of a 
massless O(N) scalar field theory with quartic interactions as given by the Lagrangian 
[58, 59, 60, 61, 62] 

£(x) = i(9^) 2 --^.g o 2 (0 2 ) 2 (3.1) 



iV + 2" 



(.2 



where there are N scalar fields fa and d> — faf + • • • 4>'n 

In the limit of TV — > 00 the only Feynman diagrams that survive are those 
that derive from ring ("daisy") diagrams or nested rings ("superdaisies" , "cactus", 
or "foam" diagrams) as shown in figure 2. 

In dimensional rcgularization, the zero-temperature scalar theory is massless 
provided the bare mass is zero. The coupling however receives an infinite 
renormalization in n — > 4 dimensions. In the large- N limit this is determined in 
modified minimal subtraction (MS) by 

113 1 

7 + ^ 



2-2+— T— • (3-2) 
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The sign of the counterterm in (3.2) in fact hints at the well-known problem of 
triviality of 4> 4 theories [61]. If the bare coupling g 2 is positive and n approaches 4 
from below, the renormalized coupling g goes to zero. As discussed in reference [62], 
if one insists on a non-trivial theory with g 2 > (which is only possible when g 2 is 
negative, and divergent for n — > 4), one finds that there is a tachyon (Landau pole) 
with mass 

m Lhyon = CX P (j£2 + 2 ) ( 3 ' 3 ) 

which appears to disqualify this model completely. Here p 2 — Aire^/j, 2 and /x is a 
mass scale introduced to make the coupling dimensionless in n =/= 4 dimensions. 

However, at small renormalized coupling, g 2 <C 1, the tachyon's mass is 
exponentially large. If everything is restricted to momentum scales smaller than (3.3), 
e.g. by a slightly smaller but still exponentially large cutoff, the n = 4 theory seems 
perfectly acceptable. For our purposes we shall just have to restrict ourselves to 
temperature scales smaller than (3.3) when considering the finite-temperature effects 
in this scalar theory. 



3.2. Thermal masses 

To one-loop order, the scalar self-energy diagram is the simple tadpole shown in the 
first diagram of figure 3, which is quadratically divergent in cutoff rcgularization, but 
strictly zero in dimensional rcgularization. The Bose distribution function occurring 
at nonzero temperature provides a cutoff at the scale of the temperature which gives 

H = (m«) 2 = 4! g 2 J ^e( q )6( q 02 q 2 ){n{f) + |} = g 2 T 2 (3.4) 

so the initially massless scalar fields acquire a temperature-dependent mass. As wc 
shall see, in more complicated theories the thermal self-energy will generally be a 
complicated function of frequencies and momenta, but the appearance of a thermal 
mass scale ~ gT is generic. 

It should be noted, however, that thermal masses are qualitatively different from 
ordinary Lorentz-invariant mass terms. In particular they do not contribute to the 
trace of the energy-momentum tensor as m 2 T 2 , as an ordinary zero-temperature mass 
would do [63]. So while the dispersion law of excitations is changed by the thermal 
medium, the theory itself retains its massless nature. 

At higher orders in perturbation theory, the thermal contributions to the scalar 
self-energy become nontrivial functions of frequency and momentum which is complex- 
valued, implying a finite but parametrically small width of thermal (quasi) particles 
[64, 65], so that the latter concept makes sense perturbatively. 



Thermal field theory 13 




Figure 3. One-loop correction to the self-energy in scalar </> 4 theory, and some- 
higher-loop diagrams with increasing degree of infrared singularity when the 
propagators arc massless 



In the large-iV limit the self-energy of the scalar field remains a momentum- 
independent real mass term also beyond one-loop order and is given by the Dyson 
equation 

n = 4! 5o V-"/ ^W^ + iJW-n). (3.5) 

(Note that in (3.4) we had simply replaced g 2 by g 2 , disregarding the difference as 
being of higher than one-loop order.) The appearance of a thermal mass introduces 
quadratic ultraviolet divergences in II, which are however of exactly the form required 
by the renormalization of the coupling according to (3.2). Including the latter, one 
finds a closed equation for the thermal mass of the form 



M / <^t»( ? °)%°)% 2 - -4) + j&i (>» > - 1) } 

=: 24 3 2 {/ T (mth) +/ / Kh,M)} • (3.6) 

The last term in the braces in (3.6), which has been neglected in [58, 66, 10], is 
responsible for a non-trivial interplay between thermal and vacuum contributions. Its 
explicit dependence on p, is such that it cancels the one implicit in g 2 — g 2 {p). 



dg 3 3 



9*, (3.7) 



d]n(fi) At: 2 ' 

which is exact for N^oo. 

More complicated effects of zero-temperature renormalization on reorganized 
thermal perturbation theories have been discussed in a scalar 0| model in [67]. Another 
solvable toy model is given by the large-iVy limit of QCD or QED, which has been 
worked out in a thermal field theory context in [68, 69, 70, 71]. These theories also 
require the introduction of a cutoff to avoid the Landau singularity and triviality. In 
contrast to the 0(iV — > oo) <ft theory, they involve complicated momentum-dependent 
dispersion laws as well as damping effects. While the large- N <j)\ theory has unphysical 
instabilities above a certain temperature (aside from the Landau singularity), large- Nf 
QED and QCD is well-behaved and very useful as a benchmark for approximations to 
real QED and QCD. 



3.3. Perturbation series 

The integral appearing in (3.6) can be evaluated e. g. using Mellin transformation 
techniques [72] to obtain a scries expansion of m t h whose first few terms read [62] 

#=/- 5 f + |;(3-7-.n i | f ) + |f(2 7 + 21» i | f - 1 ) + 0( 9 «). (3.) 

This result shows that ordinary perturbation theory is unable to go beyond the 
one- loop result (3.4), for ordinary perturbation theory is an expansion in powers of g 2 . 
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Equation (3.8) however involves odd powers of g. Indeed, because of the masslessness 
of the scalar theory, ordinary perturbation theory encounters infrared divergences 
starting at two-loop order, which are exacerbated by the Bose distribution function 
behaving as n(q°) ~ T/q° for q° — ► 0. For example, the second diagram in figure 3 
involves two massless propagators with equal momentum. The inserted tadpole does 
not vanish at zero momentum, but is given by the constant term (3.4). At zero 
temperature, such an insertion (if nonzero), would make this diagram logarithmically 
infrared divergent; at finite temperature it is instead linearly infrared divergent. The 
higher- loop diagrams shown in figure 3 are even more infrared divergent. But, as we 
have seen, the full propagator contains the thermal mass, so all these divergences are 
spurious — they just signal the need for using a resummed i.e. massive propagator. 

A systematic method to perform the required resummation of ordinary 
perturbation theory is to add a thermal mass term — l/2m 2 h cj) 2 to the Lagrangian (3.1) 
and to subtract it as a counter-term which is treated as a one-loop-order quantity. The 
corresponding calculation at N = 1, where (3.6) is only part of the full result, has 
been performed up to and including order g 5 in [64, 73] . 

3.4- Dimensional reduction 

An important technical concept for studying static quantities such as thermodynamic 
potentials and (static) screening masses is that of "dimensional reduction" [74, 75, 
76, 77, 78, 79, 80], which in the case of scalar field theories has been worked out in 
[81]. In this approach one separates hard (k <~ T) from soft (k < gT) modes and 
integrates out the former. Since in the Matsubara formalism all non-static modes 
are necessarily hard, this yields a three-dimensional effective theory containing zero- 
modes only whose parameters (masses and coupling constants) are to be determined 
by perturbative matching to the full theory. To lowest order in scalar 4 theory, this 
yields (now for N = 1) 

U = |(V0o) 2 + \m\<& + 9Ut + ■■■ (3.9) 

where </> = VT Jq dT<f>{r, x) and, to lowest order, m§ = g 2 T 2 , g 2 = g 2 T. 

Calculating now also soft one- loop corrections, one obtains in dimensional 
regularization 

«* - - -I** - <"»» 

which is exactly the correction given in (3.8). (In cutoff regularization there would be 
an extra term cx g 2 A, cancelling a contribution oc g 2 TA to m 2 .) Using this method, 
the perturbative expansion of m th for N = 1 has been worked out to order g 5 in [73] . 

3.5. Apparent convergence 

The large-TV limit provides an instructive opportunity to study the convergence 
properties of a resummed perturbation theory with the exact result obtained by simply 
solving (3.6) numerically. This has been carried out in great detail in [62] with the 
result that only for g sufficiently smaller than 1 there is quick convergence, which 
further deteriorates if the renormalization scale p, is very different from T. An optimal 
value turned out to be p, m 2ttT, the scale of the bosonic Matsubara frequencies, which 
has been argued previously in [81] to be a natural choice. 
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Figure 4. Thermal mass in l&rge-N </> 4 -thcory as a function of g(fl = 2-irT) 
together with the perturbativc results (3.8) accurate to order g 2 , g 3 , g 10 . 
The g 3 result is the one reaching zero at g 1; its Padc-improved version (3.12) 
is given by the long-dashed line. The short-dashed line just below the exact 
result is obtained by solving the quadratic gap equation (3.13), which is also 
perturbatively equivalent to the order g A result. 




However, the convergence of the resummed thermal perturbation series seems 
to be surprisingly poor given that the exact result following from (3.6) is a rather 
unspectacular function. Naturally, a perturbativc result is a truncated polynomial 
and thus bound to diverge more and more rapidly at large coupling as the order is 
increased. This may be the case even when the physical effects described by the 
lowest-order terms are still predominant. 

In [82] it has been shown that the alternative so-called nonlinear ^-expansion 
scheme [83] yields approximations to m 2 h that converge almost uniformly in g, but 
this scheme has not yet found applications in more complicated (gauge) field theories. 

A simpler alternative is provided by Pade approximants [84], that is rewriting 
a given perturbative result as a perturbatively equivalent rational function in g by 
replacing 

F n (g)=c + c l9 +... + c n g -» F [p , n _ p] (g) = ——- r ————. (3.11) 

For example, truncating (3.8) above order g 3 gives an approximation F3 to m 2 h that 
stops growing as a function of g at g — 2ir/9 ~ 0.7 and goes back to zero and then to 
negative values for g > tt/3 ~ 1 (line labelled "(3)" in figure 4). The true (large-iV) 
thermal mass from (3.6) however is a monotonically growing function (line labelled 
"exact" in figure 4). The result F3 in fact ceases to be an improvement over the 
leading-order result F2 = g 2 T 2 roughly where it stops growing. On the other hand, 
the simplest possibility for a perturbatively equivalent Pade approximant, f[2,i] ; is a 
monotonic function in g, 

»4/** = n^ + 0Grt (3.12) 

(long-dashed line in figure 4), and it gives a substantial improvement for g > 1. 

Higher-order Pade approximants converge rather well in the simple scalar model 
[62] (except when they happen to have a pole of the denominator at positive coupling) 
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and they have been proposed as a possibility to improve also the unsatisfying 
convergence of perturbative results for the thcrmodynamical potential in finite- 
temperature QCD [85, 86, 62], which will be discussed in more detail in Sect. 4.2. 
There it does increase the apparent convergence of the resummed perturbation scries 
for the first few orders, but at higher orders it looks less convincing. In fact, at these 
higher orders the perturbation scries also involves ln(g)-contributions, which make a 
simple Pade improvement appear less natural. Indeed, in the above large-TV scalar 
model, where the higher Pade approximants converge rather quickly [62], no ln(g)- 
terms arise. 

Already the simplest Pade resummation (3.12) suggests that the low quality of 
standard perturbative results is due to the fact that the latter are polynomials in g 
which inevitably blow up at larger values of g. In particular those contributions which 
can be traced to a resummation of the screening mass involve large coefficients. Since 
this resummation is a priori nonperturbative in that it involves arbitrarily high powers 
of g, this signals the need for a more complete treatment of such resummation effects. 

In the above scalar toy model, one can in fact easily obtain an efficient 
resummation of the term involving g 3, in (3.8) which is responsible for the poor 
convergence of the perturbative results displayed in figure 4. If one just retains the 
first two terms in a (m/T) expansion of the one-loop gap equation (3.6), one ends up 
with a simple quadratic gap equation 

m t \ = g 2 T 2 - -g 2 Tm th (3.13) 

7T 

which is perturbatively equivalent to (3.8) to order g 3 , but the solution of (3.13) turns 
out to be extremely close to that of the full gap equation if the MS renormalization 
scales p, w 2-kT. 

There exist formalisms which at a given order of approximations perform 
a complete propagator resummation: the so-called self-consistent ^-derivable 
approximations [87] (in particle physics also known in connection with the composite- 
operator effective action or Cornwall- Jackiw-Tomboulis formalism [88]). 

In the Luttingcr-Ward representation [89] the thermodynamic potential £1 = 
— PV is expressed as a functional of full propagators D and two-particle irreducible 
(2PI) diagrams. Considering a scalar field theory with both cubic and quartic vertices 
for the moment, Vt[D\ has the form 

Q[D] = -TlogZ = \T TrlogZ)- 1 - \T Tr HD + T$[D] 

|A n(w ) i m [logZ}- V, k) - U(u, fe)D(w, k)] + T$[£>], (3.14) 

where Tr denotes the trace in configuration space, and $[£>] is the sum of the 2- 
particle-irreducible "skeleton" diagrams 

-$[£>]= 1/120 + i/ 8 QO+i/48Q + ... (3.15) 

The self energy II = D^ 1 — Dq 1 , where Do is the bare propagator, is related to &[D] 
by 

8<f>[D]/6D = ifl. (3.16) 

An important property of the functional fl[D], which is easily verified using (3.16), is 
that it is stationary under variations of D: 

Sn[D]/5D = 0. (3.17) 
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Self-consistent ( "^-derivable" ) [87] approximations are obtained by selecting a class 
of skeletons in <&[D] and calculating II from equation (3.16) above, preserving the 
stationarity condition. 

^-derivable approximations have been worked out in scalar theory to 3-loop order 
[90, 91] and it has recently been shown that in this model these approximations can 
be nonperturbatively renormalized in a self-consistent manner [92, 93, 94, 95]. 

In the present large-TV </> 4 toy model, the 2-loop ^-derivable approximation is 
in fact exact, so that one may use this example for studying the quality of further 
approximations on top of the former [5], which are typically unavoidable to make the 
formalism tractable in more complicated theories. 

In particular, a resummation of quasi-particle propagators in thermodynamic 
quantities is conveniently performed to 2-loop order in the entropy density S — dP/dT, 
which turns out to have the simple representation [96] 



where II is the one-loop self-energy calculated self-consistently with dressed propagator 
D. At this order all fundamental interactions can be completely absorbed in the 
spectral properties of the quasi-particles, whose residual interactions enter only at 
three-loop order. For a real momentum-independent self-energy II = to 2 , (3.18) 
even coincides with the free (Stcfan-Boltzmann) expression for the entropy density, 
Sf Tee (m). The more general form (3.18) is also applicable in gauge theories and for 
fermions [97, 96, 98, 99] and we shall return to that in Sect. 6.3.2. 

In references [100, 101, 102, 103] a different reorganization of thermal perturbation 
theory has been proposed, "screened perturbation theory" (SPT), which amounts to 
adding and subtracting a mass term to the Lagrangian (3.1) according to 



and to treating <5 as a one- loop quantity prior to putting it eventually to <5 = 1. 
The difference to conventional resummation techniques is to refrain from identifying 
and expanding out the coupling g implicit in to, which in the end is chosen as some 
approximation to the thermal mass and thus proportional to g for small coupling. 

Starting from two-loop order, it is possible to determine to by a principle 
of minimal sensitivity, which makes SPT a variant of the so-called variational 
perturbation theories (see e.g. [104, 105]). 

This scheme has been applied with apparent success to scalar field theory 
[100, 101, 102, 103]. It has also been generalized to gauge theories, where a local mass 
term is insufficient but needs to be replaced by a nonlocal gauge invariant extension 
of a thermal mass term (see sect. 6.3.1) [106, 107, 108, 109]. A special difficulty of 
SPT is that at any finite order of the new perturbative expansion it gives rise to new 
ultra-violet divergences and corresponding new scheme dependences, which need to 
be fixed in some way or other [110]. 

3. 6. Restoration of spontaneously broken symmetry 

A physically important case ignored so far is spontaneous symmetry breaking, the 
simplest example of which is provided by adding a wrong-sign mass term to the </> 4 
theory considered above 




(3.18) 



C -» C - \m 2 (t> 2 + S\ 




(3.19) 



C = \{d^f - V cl (4>) with V d (<f>) = -\v 2 + 



(3.20) 
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where we have switched to the more conventional notation of A/4! = g 2 . While C is 
symmetric under <fi —> —<fi, this symmetry is "spontaneously broken" by choosing one 
of the minima of V c \, which are given by m ; n = ±y / 3\v 2 / \. 

At high temperatures, however, there is symmetry restoration [111, 112, 58, 113, 
114]: When T » the scalar field receives a contribution of rh 2 h = XT 2 / 41 to its 
(initially negative) mass squared: — v 2 — > — v 2 + AT 2 /4! or, equivalently, 

V d - V eS (T) = V d + \rh\^- (3-21) 



As a result, the minimum of V c \ becomes <pmin = for T > T c = yjA\v 2 /A. Since 

P = -K ff (T)| = ^9(T C - T) x (T c 2 - T 2 ) 2 , (3.22) 

the phase transition is of second order, i.e., there is no discontinuity in the first 
derivative of the pressure (in the entropy), but only in its second derivative, the 
specific heat. 

The above effective potential can be derived more directly and systematically from 
the partition function evaluated at spatially constant field configurations 4>, 

V eS &T) = -±lnz\_. (3.23) 

At one-loop order the temperature-dependent contribution of a scalar field with (field- 
dependent) mass m 2 {4>) = V^[{4>) is 

v£> {4>; T) = V c1 + tJ ^ In (l - e~^ T ) with ui = Vp 2 +m 2 (cf>) 

= Vcl ~^r + h m *ti )Ti nb m3( ^ )T + • • • (3 - 24) 

The first field-dependent correction term reproduces the effective potential of 
(3.21). Subsequent terms in the high-T expansion as well as higher order loop 
contributions require a resummation of the thermal mass of the scalar field: 

m 2 {$) - m c 2 ff (0; T) = m 2 (0) + m 2 h . (3.25) 

To avoid over-counting, this resummation has to take place either through the explicit 
introduction of thermal counterterms [114] or by using a self-consistent formalism like 
that of <&-derivable approximation or the CJT effective action [115]. The latter allows 
to include supcrdaisy diagrams, which treated naively would give a spurious and in 
fact completely misleading contribution to the effective potential oc m((p)T 3 . 

Insertion of (3.25) in (3.24) gives an effective potential which at the critical 
temperature has two degenerate minima, = and cfr ~ \/AT c . However, at the 
second non-trivial minimum, m c g(<fi) ~ AT and the loop expansion parameter of the 
perturbation series, which is XT/m e s, ceases to be small. Hence, perturbation theory 
cannot decide whether there is a first-order phase transition as implied by a second 
degenerate minimum in the potential or not. In fact, universality arguments show 
that the phase transition must be second order for a </> 4 theory with Z 2 symmetry 
[61]. 

The situation is different in gauge theories. Again, resummed perturbation 
theory signals a first order phase transition, but when the Higgs mass is small 
compared to the massive vector-boson mass, it is sufficiently strongly first-order so 
that the perturbation series does not break down where the degenerate minima appear. 
The (resummed) effective potential has been calculated pcrturbatively to one-loop 
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[116, 117, 118] and two-loop order [114, 119, 120] (see also [121]). Of particular interest 
is the question whether the electroweak sector of the Standard Model (or extensions 
thereof) admits a first order phase transition, which would be of great interest for 
baryogenesis in the Early Universe [122]. Nonperturbative studies indeed confirmed 
a first-order transition for sufficiently small Higgs masses [123, 124], but found an 
cndpoint [125, 126, 127] at Standard Model Higgs mass ran < 80 GeV, so that a first- 
order phase transition in the Standard Model is excluded by the current experimental 
bounds on the Higgs mass (though not in possible (supersymmetric) extensions of the 
Standard Model [128, 129]). 

4. QCD thermodynamics 

The (resummcd) perturbative evaluation of the thermodynamic potential of QCD 
at high temperature has been pushed in recent years up to the order <? 6 In g 
[130, 131, 132, 133, 134, 135, 136, 137]. At higher orders in g this is much facilitated 
by the possibility to employ effective field theory methods which in this case lead 
to a dimensional reduction to a 3-dimcnsional Yang-Mills plus adjoint Higgs theory 
[75] . Completion of these results at order g 6 is in fact impossible without inherently 
nonperturbative input, but further progress has been made most recently by the 
extension to nonzero quark chemical potential within dimensional reduction [138]. 

At zero temperature and high chemical potential the pressure is known to order 
g A [139, 140, 141, 142, 138], while the low-temperature expansion of the pressure leads 
to the phenomenon of non- Fermi-liquid behaviour of entropy and specific heat [143]. 

4-1. Dimensional reduction 

Dimensional reduction in hot QCD leads to an effective three-dimensional Lagrangian 
[75, 76, 77] 

C E = \ tr + tr[A, A)] 2 + m\ tr A\ + ±A B (tr A\f + . . . (4.1) 

where the parameters are determined perturbatively by matching [136, 144]. In lowest 
orderS and at zero chemical potential one has a dimensionful coupling g\ = g 2 T and 
[77] 

9 — Nc 

m| = (l + A7/6). 9 2 T 2 , X E = —^g*T, (4.2) 

1Z7T 

though \e starts to contribute to the pressure only at order g 6 . At this order, however, 
the self-interactions of the massless magnetostatic gluons start to contribute, and these 
contributions are inherently non-perturbative because the three-dimensional theory for 
the zero modes Ai(x) is a confining theory [149, 150, 151]. 

The thermal pressure of the 4-dimensional theory can be decomposed into 
contributions from the hard modes ~ T, calculable by standard perturbation theory, 
and soft contributions governed by (4.1) which involves both perturbatively calculable 
contributions up to order g 5 T 4 and the nonperturbative ones starting at order g 6 T 4 . 

In reference [136] the effective theory based on (4.1) has been used to organize 
and reproduce the perturbative calculation of the thermal pressure to order g 5 of 
references [133, 134]. This turns out to be particularly elegant when dimensional 

E Some higher-dimension terms in the effective theory (4.1) have been determined in references [78, 
145, 146], and, including the electroweak sector, in references [147, 148]. 
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regularization is used to provide both the UV and IR cutoffs of the original and 
effective field theories. 

To order g , the contribution of the hard modes can then be written as [136] 



871^ rpj 

45 1 



(i + ¥M f (i + 



+ |244.9 + 17.247V/ - 0.415iV) + 135 (l + ±N f ) In 
-Jf(H-*N,)(l-iW/)lni}(i)*}. 



2ttT 



(4.3) 

In the first logarithm the dimensional regularization scale £l is associated with 
regularization in the infrared and thus has to match a similar logarithm in the effective 
theory, whereas the second logarithm is from UV and involves the first coefficient of 
the beta function. 

Indeed, calculating the pressure contribution of the soft sector described by (4.1) 
in dimensional regularization gives, to three-loop order (neglecting A^-contributions) 



PmA /T = A m | - ^ 
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(4.4) 



All the contributions to the pressure involving odd powers of g in (4.5) (as well as 
part of those involving even powers) are coming from the soft sector. Inserting the 
leading-order value (4.2) for ttle gives the QCD pressure up to and including order 
g 4 In 5; to obtain all the terms to order g 5 , ncxt-to-leading order corrections to the mj- 
parameter have to be obtained by a matching calculation as given in reference [136]. 
The result is known in closed form [133, 134, 136], but we shall quote here only the 
case of SU(3) with Nf quark flavours and numerical values for the various coefficients: 
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(^) 5/ \o(al\na s )\. (4.5) 



Here fl is the renormalization scale parameter of the MS scheme and a s (p.) is the 
corresponding running coupling. 

The coefficient of the In a s term, the last in the pressure at high T and 
vanishing chemical potential that is calculable completely within perturbation theory, 
has recently been determined as [137, 152] 



g 6 lnfl 



8tt^ rp4 
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(4.7) 



In order to obtain the g 6 In g-contribution (4.6) one also needs g 2 E to order g 4 
(given in reference [144]) and above all the four- loop contribution of the effective 
theory (4.1) which has recently been calculated analytically as [137] 

N > { W [(« - + (* - >»4 + 

up to a constant c that is strictly nonperturbative and needs to be determined 
by three-dimensional lattice calculations. Such calculations have been undertaken 
in reference [153], but they depend on an as yet undetermined 4-loop matching 
coefficient. At the moment the conclusion is that it is at least not excluded that 
the lattice results based on dimensional reduction can be matched to the full four- 
dimensional results at temperatures of a few times the transition temperature. For 
this reason, the most reliable results on the thermodynamics of hot QCD (particularly 
for pure-glue QCD) remain to date the four-dimensional lattice data. However, since 
inclusion of fermions is particularly easy in the dimensional reduction method, but 
computationally expensive in lattice gauge theory, a full three-dimensional prediction 
would clearly be most desirable. 



4-2. Apparent convergence 
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Figure 5. Strictly pcrturbative results for the thermal pressure of pure glue QCD 
normalized to the ideal-gas value as a function of a s (p, = 2irT). 

Figure 5 shows the outcome of evaluating the perturbative result (4.5) for the 
thermodynamic pressure at Nf = to order a s , o^J 2 , a,, and a J 2 ', respectively, 
with a choice of p, — 2ttT. Apparently, there is no convergence for a s > 0.05 which 
in QCD corresponds to T < 10 5 T C , where T c is the deconfinement phase transition 
temperature. What is more, the numerical dependence on the renormalization scale 
j2 does not diminish as the order of the perturbative result is increased, but becomes 
more and more severe, as shown in figure 6, where a s (p.) is determined by a 2- 
loop renormalization group equation and p, is varied between irT and AttT. So 
there seems to be a complete loss of predictive power at any temperature of interest 
[133, 134, 136, 135]. 

To alleviate this situation, various mathematical extrapolation techniques have 
been tried, such as Pade approximants [85, 86, 156], self-similar approximants [157], 
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Figure 6. Strictly pcrturbative results for the thermal pressure of pure glue 
QCD as a function of T/T c (assuming T c /Aj^g = 1.14). The various gray bands 
bounded by differently dashed lines show the pcrturbative results to order g 2 , 
<; 3 , g 4 , and g 5 , using a 2- loop running coupling with MS renormalization point 
p, varied between ttT and 47rT. The thick dark-grey line shows the continuum- 
extrapolated lattice results from reference [154]; the lighter one behind that of a 
lattice calculation using an RG-improved action [155]. 

and Borel resummation [158, 159], however with limited success. While in the scalar 
toy model of Sect. 3.5, Pade approximants work remarkably well, in the QCD case 
there is a problem how to handle logarithms of the coupling and, perhaps related to 
that, the numerical results obtained so far appear less satisfactory. 

When compared with the perturbative results, it is however remarkable that the 
next-to-leading result to order g 2 performs rather well at temperatures > 2T C , though 
the higher-order results prove that perturbation theory is inconclusive. Moreover, 
simple quasiparticle models which describe the effective gluonic degrees of freedom by 
2N g (Ng = N 2 — 1) scalar degrees of freedom with asymptotic masses taken from a 
HTL approximation can be used quite successfully to model the lattice data by fitting 
the running coupling [160, 161, 162, 163]. 

In fact, we have seen in section 3.5 that already in the simplest scalar model 
resummed perturbation theory gives rather poorly convergent results, and that fairly 
simple reorganizations as in (3.13) lead to dramatic improvements of the situation. 

As has been noted in references [137, 164], in the calculation as organized through 
the dimensionally reduced effective theory (4.1), the large scale dependence of strict 
perturbation theory can be significantly reduced when the perturbative values of 
the effective parameters are kept as they appear in (4.4), without expanding and 
truncating the final result to the considered order in the coupling. What is more, 
the unexpanded two-loop and three-loop contributions from the soft sector lead to 
results which no longer exceed the ideal-gas result (as the strictly perturbative results 
to order g 3 and g A do), and their (sizable) scale dependence diminishes by going from 
two to three-loop order [164]. 

At three-loop order, it is in fact possible to eliminate the scale dependence 
altogether by a principle of minimal sensitivity. The result in fact agrees remarkably 
well with the 4-d lattice results down to ~ 2.5T C as shown in figure 7. 

The four-loop order result depends on the unknown constant c in (4.7). However, 
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Figure 7. Three-loop pressure in pure-glue QCD with unexpanded effective-field- 
theory parameters when p is varied between izT and 4ttT (medium-gray band) ; the 
dotted lines indicate the position of this band when only the leading-order result 
for niE is used. The broad light-gray band underneath is the strictly perturbative 
result to order g 5 with the same scale variations. The full line gives the result 
upon extremalization (PMS) with respect to p (which does not have solutions 
below ~ 1.3T C ); the dash-dotted line corresponds to fastest apparent convergence 
(FAC) in m|, which sets p f» 1.79irT. (Taken from [164]) 
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Figure 8. Like figure 7, but extended to four-loop order by including the recently 
determined </ 6 ln(l/</) contribution of [137] together with three values for the 
undetermined constant 8 in [g 6 ln(l/cj)+<5]. The broad light-gray band underneath 
is the strictly perturbative result to order g 6 corresponding to the central value 
<5 = 1/3, which has a larger scale dependence than the order g 5 result in figure 
7; the untruncatcd results on the other hand show rather small scale dependence. 
The full line gives the untruncated result with 5 = 1/3 and p fixed by PMS 
(which does not have solutions below ~ 1.9T C ); the dash-dotted line corresponds 
to fastest apparent convergence (FAC) in m 2 E , which sets again to p s» 1.79-7rT. 
(Taken from [164]) 




it is at least not excluded that c could be such that the four-loop result is also close to 
the 4-d lattice results [137, 165]. In fact, while a strictly perturbative treatment leads 
to increased scale dependence compared to 3-loop order, keeping the soft contributions 
unexpanded in g further diminishes the scale dependence [164] as shown in figure 8. 
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While this goes only minimally beyond a strictly perturbative treatment, it 
strongly suggests that perturbative QCD at high temperature, when supplemented 
by appropriate resummation of soft physics, is not limited toT> 10 5 T C as previously 
thought [133, 135], but seems capable of quantitative predictions at temperatures of 
possibly only a few times the transition temperature. 



4-3. Finite chemical potential 

The calculation of the thermodynamical potential using dimensional reduction can in 
fact be extended to nonvanishing quark chemical potential /j, q , provided however that 
T 3> rriE- This has been carried out up to and including order g 6 In g in reference [138], 
and a discussion of the issue of apparent convergence can be found in reference [166]. 

A non-zero chemical potential modifies the parameters of the effective theory, in 
particular the mass parameter to_e, which is altered already at leading order according 
to (5.18) below. In addition there are new, C-odd terms in the effective Lagrangian. 
The one with smallest dimension in nonabelian theories reads [167, 168, 169] 

<? 

In Abelian theories there is also a linear term involving tr A which has been discussed 
e.g. in [170]. 

In general the effects of these additional C-odd terms are small compared to 
the C-even operators in (4.1) which depend on the chemical potential through its 
parameters. One quantity which is determined to leading order by the operator (4.8) 
is the flavour off-diagonal quark number susceptibility at zero chemical potential [171] 

d 2 P 

Xij = t-JT- ■ ( 4J ) 

When quark masses are negligible, all off-diagonal components are equal at /Xj = 0. 
Denoting them by x, the leading-order term involves a logarithmic term coming from 
the exchange of three electrostatic gluons and is given by [171] 

(N 2 - 1)(N 2 - 4) /g\6 „ 1 

where N is the number of colours. This vanishes in SU(2) gauge theory, but not in 
QED, where (in the ultrarelativistic limit) [171] 

„6 i 

X - — *T 2 \n-. (4.11) 

A qed 24tt 6 e V ' 

On the other hand, the diagonal quark susceptibilities have a perturbative 
expansion whose first few terms are given by 



X i3N g(g y 3N n N f (g y 3N g , g y ,1 , ofo4l (4U] 

with xo = NT 2 / '3 the ideal gas value and N g = N 2 — 1. The higher-order coefficients 
have been calculated by Vuorinen [172] up to and including order g 6 \n(l/g). The 
problem with apparent convergence is similar if somewhat less severe than in the case 
of the pressure at zero chemical potential discussed above. 

In contrast to the pressure, however, the coefficient of the order g 6 term 
in x is n °t sensitive to nonperturbative chromomagnetostatic physics and thus 
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calculable in perturbation theory, though not yet available. Its determination would 
in fact be of some interest in view of the important progress that has recently 
been made with the inclusion of small chemical potentials in lattice gauge theory 
[173, 174, 175, 176, 177, 178, 179, 180]. 



4-4- Low temperature and high chemical potential 

At large chemical potential but T < g/j, q , dimensional reduction does not occur. In 
this case a fully four-dimensional computation has to be performed. 

The perturbative result up to and including order g A at zero temperature has 
been calculated by Freedman and McLerran [139, 140, 141] and by Baluni [142] more 
than a quarter of a century ago. 

This result was originally given in a particular gauge-dependent momentum 
subtraction scheme. In order to convert it to the gauge-independent MS scheme, one 
needs to replace the scale parameter fio in reference [141] (M in references [142, 8]) 
according to [99, 181] 

=Mexp{[(151)7V-407V/]/[24(117V-27V / )]}. (4.13) 

Furthermore, the order g A contributions involved two integrals that were evaluated 
only numerically with sizable error bars. This calculation was recently repeated in 
reference [138], with one of these integrals evaluated analytically and the other with 
very high accuracy. 

Specialized to N — 3 and uniform chemical potential for Nf quark flavours, the 
result for the pressure at zero temperature to order g 4 in the MS scheme reads 



"—iP* 1 - 3 - 



18 - 11 log 2 - 0.53583AT/ + N f log ■ 



+ (n-lN f )lo g ^}^r + 0(al\o g a s )\, (4.14) 

where fx q is the (common) quark chemical potential, not to be confused with the 
renormalization scale p, of the MS scheme. 

At zero temperature and large chemical potential, perturbation theory is not 
hampered by a perturbative barrier at order g 6 so that higher-order corrections are 
in principle calculable, but not yet available (except for one term oc Njg 6 extracted 
numerically from the solvable large- Nf limit of QCD [71]). 

However, at small but finite temperature T < g/j,, the only weakly screened low- 
frequency transverse gauge-boson interactions (see (5.19) below) lead to a qualitative 
deviation from the Fermi liquid behaviour of relativistic systems described in reference 
[182]. In particular, the low-temperature limit of entropy and specific heat does not 
vanish linearly with temperature, but there is a positive contribution proportional 
to aT In T _1 [183, 184, 185], which implies that at sufficiently small temperature 
the entropy exceeds the ideal-gas result. In reference [143] this effect was most 
recently calculated beyond the coefficient of the leading log obtained in [183, 185] ||. 
The complete result for the entropy below order T 3 In T, where regular Fermi- liquid 
corrections enter, is given by (for SU(3) and with numerically evaluated coefficients 

|| The coefficient of the aTlnT -1 term given in the original paper [183] was found to be lacking of 
a factor of 4. 



Thermal field theory 



26 



which arc known in closed form [143]) 



I l 

+ 0.130147V S f^-V f— Vl+0(T 3 lnT), T « g N , (4.15) 



and turns out to involve also fractional powers of the temperature as well as the 
coupling. (The corresponding result in QED is obtained by replacing N g —* 1 and 
a s Nf — > 2a.) Equation (4.15) is the beginning of a perturbative expansion provided 
T/n q <C g, e.g. T//i q ~ g 1+s with <5 > 0. The corresponding contribution to the 
pressure is then of the order g 4+2S In g and thus of higher order than the terms 
evaluated in (4.14). But in the entropy the zero- temperature limit of the pressure drops 
out and the nonanalytic terms found in (4.15) become the leading interaction effects 
when T I [i q <C g. For exponentially small T/{g^L q ) <~ cxp(— #/.g 2 ), they eventually 
become comparable to the ideal-gas part and the perturbative treatment breaks down 
(in QCD one in fact expects a breakdown of perturbation theory already at the order 
of exp(—#/g) for colour superconducting quarks, cf. section 5.3). 

Originally the anomaly in the specific heat was discussed for a nonrelativistic 
electron gas with the expectation that this effect may be too small for experimental 
detection [183]. In QCD, however, it is numerically much more important, not 
only because a s » QqeDj but also because of the relatively large factor N g = 8. 
Consequently, it may play some role in the thermodynamics of (proto-)neutron stars, 
if those have a normal (non-superconducting) quark matter component. 

For potential phenomenological applications in astrophysical systems, the specific 
heat C v at constant volume and number density is of more direct interest, which 
however differs from the logarithmic derivative of the entropy only by sublcading 
terms [186]: 

where N is the number density. In QCD this gives 



Cv = + ~"i8n y ln 1 L9574 \/^^ ) - 0.288095iV 3 



JMV> , a s NfN g ^ n _ , la s N f N \ n oaatMtlKr fasN£\ * ( T_ 

fJ-q 



+ 0.3036697 N g (^^y (Jpj* J +(9(T 3 lnT), T < gfiq . (4.17) 

The results (4.15) and (4.17) imply an excess over their respective ideal-gas values. 
They depend on having T -c g^ q , and in this regime the low-temperature series has 
a small expansion parameter T/(g/j, q ); for T ^> g\i q the standard exchange term gives 
the leading-order interaction contribution [8] 

Cv* S = ^ q T |^ - N g N f %- + 0(ai)} + 0(T 3 ), (4.18) 

so that for larger temperature there is a reduction compared to the ideal-gas result. 
A nonmonotonic behaviour of the entropy as a function of T which interpolates 
between (4.15) and (4.18) has indeed been found in the numerical evaluation of the 
exactly solvable large- Nf limit of QED and QCD [70], and there the domain where 
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the entropy has the anomalous feature of exceeding the ideal-gas value is given by 



Non-Fermi-liquid corrections to C v in the context of ultrarelativistic QED and 
QCD have also been considered in reference [187], however the result obtained therein 
does not agree with (4.17). When expanding pcrturbatively the renormalization-group 
resummed result of reference [187] , it would imply a leading nonanalytic aT 3 In T term, 
which is in fact the kind of nonanalytic terms that appear also in regular Fermi-liquids 
[188]. However, reference [187] did not evaluate all contributions oc aT, which were 
considered to be free of nonanalytic terms. 

5. The quasi-particle spectrum in gauge theories 

With the exception of static quantities in the high-temperature limit, where 
dimensional reduction is applicable, a systematic calculation of observables in thermal 
field theory, such as reaction rates, transport coefficients, and even the thermodynamic 
potential at low temperatures and high chemical potential require the determination 
and consistent inclusion of medium effects on the dynamical propagators in the theory. 

As we have seen in the example of scalar field theory in Sect. 3.2, the interactions 
with the particles of the heat bath modify the spectrum of elementary excitations. 
The poles of the propagator, which determines the linear response of the system 
under small disturbances, receive thermal corrections which, in general, introduce a 
mass gap for propagating modes and screening for non-propagating ones, even when 
the underlying field theory is massless. The simple 0(N — * oo) </> 4 model considered 
in Sect. 3.2 is in fact somewhat misleading in that there the self-energy is a real 
quantity to all orders in the coupling, whereas in a nontrivial quantum field theory 
a nonvanishing width for propagating modes is unavoidable [189]. Nonetheless, the 
typical situation in perturbation theory is that the width is parametrically smaller than 
the thermal mass at a given momentum, and that the former can be treated to some 
extent perturbatively. It should be noted, however, that such quasiparticle excitations 
need not correspond to simple poles on the unphysical sheet. There could instead be 
branch points, branch singularities, essential singularities, or even no singularities at 



Moreover, in the case of gauge theories it is a priori not clear whether the 
thermal corrections to the various propagators encode physical information or not. 
In Abclian gauge theory, the photon propagator is linearly related to correlators 
of the gauge-invariant electromagnetic field strength and so has a direct physical 
interpretation; matter fields on the other hand already transform non-trivially 
under gauge transformation and indeed their propagator is a gauge-fixing dependent 
quantity. In non-Abelian gauge theories, the gauge bosons carry colour charge and 
their propagator is also gauge dependent. Correlators of the non-Abelian field strength 
are not gauge independent cither. 

Still, even in gauge-dependent quantities there may be gauge-independent 
information. Indeed, in [191, 192] (see also [193] for a more detailed recent review) 
it has been shown that the singularity structure (location of poles and branch 
singularities) of certain components of gauge and matter propagators are gauge 
independent when all contributions to a given order of a systematic expansion scheme 
are taken into account. 

Another example for gauge-independent content in gauge-dependent quantities is 
provided by the high-temperature limit of self-energies and more-point correlation 




'/30. 



all [190]. 
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functions with small external momenta. Those are related to forward-scattering 
amplitudes [194, 195] of on-shell plasma constituents and are therefore completely 
gauge independent. They can form the building blocks of an effective theory at soft 
scales (with respect to the temperature), as we shall describe further below. Before 
doing so, we review the structure of the various propagators in a gauge theory. 



5.1. Gauge-boson propagator 

The self-energy of gauge bosons is a symmetric tensor, the so-called polarization tensor, 
which is defined by 

-lP' y = G~ lAiiy - Gq 1 ^. (5.1) 

(In the real-time formalism, where this should be defined first as a 2 x 2 matrix relation, 
we assume that II has been extracted after diagonalization, for example in the F, F 
basis according to (2.27).) 

If the gauge fixing procedure does not break rotational invariance in the plasma 
rest frame, this Lorentz tensor can be decomposed in terms of four independent tensors 
and associated structure functions, 

-ip ,y = ii a a^ + n B B» v + n c c^ + n D D» v , (5.2) 

with, in momentum space, being the spatially transverse tensor (2.33), and the 
others chosen as [131, 151, 196, 197, 11, 55] 

ft^ft 1 ' k^k v 

B» v {k) = iLJ- = rT ~ ~ A^(k), (5.3) 
n z k z 

C^ v (k) = ^-{n^ + k^n"}, (5.4) 

k^k v 

D ^( k ) = — ( 5 - 5 ) 

where n M = (r)^ a — k^k a /k 2 )5® in the plasma rest frame. 

Only A and B are transverse with respect to the four-momentum k, A^ v k v = 
B^ v k v = 0. C obeys the weaker relation C^k^k„ — 0, and D projects onto k. 
This particular basis is a convenient choice because A, B, and D are idempotent and 
mutually orthogonal; C is only orthogonal to A, but its product with the other tensors 
has vanishing trace. 

In Abclian gauge theory with linear gauge fixing, the Ward identities imply 
transversality of the polarization tensor, W^ v k u = 0. In non- Abclian gauge theories 
transversality holds only in certain gauges such as axial gauges [198] (with only 
the temporal axial gauge respecting rotational invariance) and background-covariant 
gauges [199, 200, 201]. Contrary to the experience at zero temperature, at finite 
temperature the polarization tensor turns out to be non-transverse already at one- 
loop order in general covariant and Coulomb gauges [197, 202], with the fortuitous 
exception of Feynman gauge (at one-loop order) . 

For a gauge breaking Lagrangian (2.34) BRS invariance only requires that 
[203, 11, 204, 205] 

fJuG^ = -a . (5.6) 

This entails that [204, 205] 

n D (k 2 - U B ) = n 2 c (5.7) 
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so that generally the polarization tensor contains 3 independent structure functions. 
Equation (5.7) also implies that at one-loop order II D = 0, but not beyond in those 
gauges where n c ^ 0. 

Rotationally invariant gauge fixing vectors can be written generally as / M = 
f3(k)k ft +7(/c)n A ' with (3 ^ 0. This includes covariant gauges {$ = 1,7 = 0), Coulomb 
gauges (f3 — n 2 ,7 = —A; ), or temporal gauges (/3 = /c°/fc 2 ,7 = 1). 

For these, the structure functions in the full propagator 

-G^ = A A A" V + AbB"" + A C C^ + A D D» V (5.8) 
are determined by 

A A = [k 2 -n A }- 1 (5.9) 

(3*k 2 -aIl D V ' 

5.1.1. Gauge independence of singularities In [191, 192] it has been shown that 
under variations of the gauge fixing parameters (in our case a, (3, 7) one has "gauge 
dependence identities" (generalized Nielsen identities [206, 207, 208, 209]) which are 
of the form 

SA^(k) = A^ 1 \-AZ(k)6X"(kj\ = A^(k)SY(k), (5.13) 



SA B \k) = A 



B 



n" 7/?-an c /|fc| fcM 
n 2 (3 2 k 2 - aU D 



2n„6X^ = A B 1 {k)5Z{k), (5.14) 



where 5X" has a diagrammatic expansion which is one-particle-irreducible except for 
at most one Faddeev-Popov ghost line. No such relation exists for Ac or A^. 

Now if 5Y and SZ are regular on the two "mass- shells" defined by A A X = 
and A^ 1 = 0, the relations (5.13,5.14) imply that the locations of these particular 
singularities of the gluon propagator are gauge fixing independent, for if A^ 1 = = 
A^ 1 then also A^ 1 + SA^ = = A^ 1 + 5A B \ 

In the case of As, singularities in SZ(k) include a kinematical pole 1/k 2 hidden in 
the n's and the manifestly gauge-dependent Ajj (cf. (5.12)). Excluding these obvious 
gauge artefacts, everything depends on whether the possible singularities of 6X" could 
coincide with the expectedly physical dispersion laws A^ 1 = and A^ 1 = 0. Because 
5X" is one-particle reducible with respect to Faddeev-Popov ghosts, the singularities 
of the latter have to be excluded, too. However, these are generically different from 
those that define the spatially transverse and longitudinal gauge-boson quasi-particles. 
Indeed, in leading-order thermal perturbation theory the Faddeev-Popov ghost self 
energy and the physical self energies receive contributions carrying different powers 
of temperature or chemical potential and therefore have independent and generically 
different dispersion laws. 

However, SX" may develop singularities also from onc-particle-irreducible 
subdiagrams, namely when one line of such a diagram is of the same type as the 
external one and the remaining ones are massless. This can give rise to infrared 
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Figure 9. Location of poles in and Ag of the hard-thcrmal-loop gauge 
propagator. The right part with k 2 > corresponds to propagating normal 
modes, the left part to (dynamical) screening. 



or mass-shell singularities and seemingly constitute an obstruction to the gauge- 
independence proof [210, 211]. But such singularities will be absent as soon as an 
overall infrared cut-off is introduced, for example by restricting everything to a finite 
volume first. In every finite volume, this obstruction to the gauge- independence proof 
is then avoided, and A^ 1 = and A^ 1 = define gauge-independent dispersion laws 
if the infinite- volume limit is taken last of all [212]. 



5.1.2. Gauge- field quasi-particles at leading order It turns out that the leading-order 
contributions to LT 11 ' for small frequencies and momenta k , |k| -C T are entirely 
transverse (lie = IId =0) and gauge-parameter independent. They are generated 
by one-loop diagrams with hard loop momentum and for this reason they are termed 
"hard-thermal-loops" (HTL) [213]. Their form is the same in Abclian [214, 215] and 
non-Abelian gauge theories [216, 217], and also in the presence of chemical potentials 
fif for fermions such that max(T, /x/) » k°, |k|. This universal result reads 

ii A = i(n/-n B ) (5.15) 

fi B = -^n 00 (5.16) 

with 

V = *£. ft oo = ^( 1 -^ ln F^)> ( 5 - 1? ) 

2 f e ¥/3 + £ W for QED, 

m ° ~ \ g 2 NT 2 /3 + E/5 2 M//( 2 ^ 2 ) for SU(iV) with N f flavours. ^' 18J 

The spectral representation of the resulting propagators Aa,b involves simple poles 
and continuous parts for k 2 < 0, which are given in detail in appendix A.l. 

In figure 9 the location of the poles of the corresponding propagators Aa and 
Ab are displayed in quadratic scales. The light-cone is marked by a dashed line. A 
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simple mass hyperboloid would be given by a line parallel to the latter. Evidently, 
the dispersion laws of the HTL quasiparticle excitations are more complicated — they 
involve momentum-dependent masses: in the long-wavelength limit |k| — > there is a 
common lowest (plasma) frequency u> p i. = tho/VS for propagating normal modes. 

For larger frequencies and momenta ui, |k| 3> rhp it turns out that mode A 
approaches asymptotically a mass hyperboloid with mass = mo/V^. So in this 
momentum region the physical spatially transverse polarizations of the gauge bosons 
acquire indeed a constant thermal mass. The additional mode B, whose dispersion 
curve approaches the light-cone exponentially, is found to have a spectral strength 
(the residue of the corresponding pole) that decays exponentially as k/rho — * oo 
[218], showing its exclusively collective nature. 

For real |k| but u 2 < k 2 , LI^ has an imaginary part ~ rh 2 D from the logarithm 
in (5.17) which prevents the appearance of poles in this region. This imaginary part 
corresponds to the possibility of Landau damping, which is the transfer of energy from 
soft fields to hard plasma constituents moving in phase with the field [219, 6] and is 
an important part of the spectral density of HTL propagators. At higher, subleading 
orders of perturbation theory, it is, however, not protected against gauge dependences 
in nonabelian gauge theories, in contrast to the location of the singularities which 
determine the dispersion laws of quasi-particles. However, at asymptotically large 
times, Landau damping is (generically) dominated by the gauge-independent location 
of the branch cuts at ui = ±|k|, resulting in power-law relaxation of perturbations 
[220]. There is also an exponential component of Landau damping due to a pole 
at purely imaginary u> and real |k|, on the unphysical sheet reached by continuation 
through the branch cut between w = ±|k| [221]. 

For real w < <2> p i., there are no poles for real |k|, but instead for imaginary 
Vk 2 = in, corresponding to exponential (dynamical) screening of (time-dependent) 
external sources. These poles, displayed on the left part of figure 9, are in fact closely 
related to the just mentioned poles at purely imaginary uj and real |k|. 

In the static limit lo — > 0, only mode B is screened with (Debye) screening length 
m^ 1 . This corresponds to exponential screening of (chromo-) electrostatic fields. The 
transverse mode A on the other hand is only weakly screened when ui <C mo with a 
frequency-dependent inverse screening length [217] 



When u — > 0, mode A describes magnetostatic fields which are found to be 
completely unscreened in the HTL approximation. In QED the absence of magnetic 
screening is intuitively clear and can be proved rigorously to all orders of perturbation 
theory [215, 222], but not in the nonabelian case. 

In QCD, the absence of a magnetostatic screening mass causes problems for 
perturbation theory as will be discussed further in Sect. 7.4. In fact, lattice simulations 
of gauge fixed propagators in nonabelian theories do find a screening behaviour in the 
transverse sector, though the corresponding singularity is evidently quite different 
from a simple pole [223]. 

5.2. Fermions 




(5.19) 



The fermion self-energy at non-zero temperature or density has one more structure 
function than usually. In the ultrarelativistic limit where masses can be neglected, it 



Thermal field theory 32 

can be parametrized by 

E(w, k) = a(u, |k|) 7 ° + b(u, |k|)k • 7. (5.20) 

(For a massive fermion, this would also include a mass correction, i.e., E = 07° + bk ■ 
7 + cl.) This can be rewritten as: 

7o S(^,k) = S+(w,|k|)A+(k) - E_(w,|k|)A_(k), (5.21) 

where E± = b ± a, and the spin matrices 

1 ± • k 

A ± (k) = (5.22) 

project onto spinors whose chirality is equal (A + ), or opposite (A_), to their helicity. 
Dyson's equation S^ 1 = — ft + E then implies 

7o5 ,_1 = A^ 1 A + + Al 1 A_, (5.23) 
with Aljj 1 = — [o;=F(|k| + E±)]. This is trivially inverted to yield the fermion propagator 

S70 = A+A+ + A_A_. (5.24) 
whose singularities are conveniently summarized by the equation 

detS _ V,k) = 0, (5.25) 

where the determinant refers to spinor indices. The potential gauge dependences of 
the singularities of S are described by an identity of the form [192] 

<5dctS -1 (o;,k) =detS -1 (o;,k) <5trX(w,k), (5.26) 

where 5tr X(u;,k) again has a diagrammatic expansion which is one-particle- 
irrcducible except for at most one Faddeev-Popov ghost line. The same reasoning 
as in Sect. 5.1.1 (with similar qualifications) leads to the conclusion that the positions 
of the singularities of the fermion propagator are gauge- fixing independent [192]. 

In the HTL approximation, the fermion self-energies are once again gauge- 
independent in their entirety. Explicitly, they read [224, 225, 226, 227, 228, 229]: 

where M 2 is the plasma frequency for fermions, i.e., the frequency of long- wavelength 
(k — > 0) fermionic excitations: 

M 2 = t^L (t 2 + 4 V (5.28) 



8 

(C f = (N 2 - 1)/2N in SU(A) gauge theory, and g 2 C f -» e 2 in QED.) 

For frequencies lo < M, there arc, in contrast to the gauge boson propagator, no 
solutions with imaginary wave-vectors that would correspond to screening. Instead, 
the additional collective (— ) (occasionally dubbed "plasmino" [230]) branch exhibits 
propagating modes down to w 0.928M (at |k|di p ~ 0.408M), with a curious dip 

in the dispersion curve reminiscent of that of rotons in liquid helium [231]. 

For momenta |k| 3> M, the normal (+) branch of the poles of the fermion 
propagator approaches asymptotically a mass hyperboloid with mass 
and unit residue, whereas the (— ) branch tends to the light-cone exponentially, with 
exponentially vanishing residue. 




For space- like momenta, ui 2 < k 2 , there is again a large imaginary part oc M 2 
corresponding to fcrmionic Landau damping, which now reflects the possibility of 
turning hard fcrmionic (bosonic) plasma constituents into hard bosonic (fermionic) 
ones in the presence of soft fermionic fields. Beyond lowest (HTL) order, this imaginary 
part cannot be expected to be gauge-fixing independent, however, not even in Abelian 
theories. 

Another case where gauge dependences arise is when the bare fermion mass cannot 
be neglected compared to temperature. Then the HTL approximation is no longer 
adequate, and one has to include more than one-loop diagrams for consistency. At 
one-loop order, the modifications of the fermionic dispersion laws have been studied 
in references [232, 233], with the result that the additional (— ) branch disappears 
gradually when m/T > 1/3. This phenomenon has been obtained both in Coulomb 
and Feynman gauge, with only weak gauge dependences in the real parts of the fermion 
self-energy, so it is not to be expected that higher-loop corrections would change this 
result qualitatively. 

There have been also investigations of how the HTL dispersion laws of fermions 
as well as gauge bosons are affected by retaining non-leading powers of temperature 
[234, 235] and chemical potential [236, 237, 238, 239] which also emphasize the need 
for higher-loop contributions to obtain complete, gauge-independent results. 

The very existence of the additional (— ) branch in the high-temperature limit 
has been confirmed in references [240, 241] under rather weak assumption without use 
of perturbation theory. Possible experimental signatures of the additional fcrmionic 
quasi-particles are the Van Hove singularities in dilepton production from a quark- 
gluon plasma [242, 243, 240] (though smeared out by damping effects to an unknown 
extent); they also play a role in the calculation of the electroweak baryon genesis of 
reference [244] but it has been shown [245, 246] that the fermion damping discussed 
in section 7.5 completely swamps their effects in the standard model (but may be still 
important in extensions thereof [247]). 
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5.3. Diquark condensates and colour superconductivity 

At sufficiently low temperature and high quark chemical potential, QCD will in fact 
be in a colour superconducting phase which modifies the fcrmion propagator by the 
appearance of a diquark condensate. The reason is that even an arbitrarily weak 
attractive interaction at the Fermi surface leads to the appearance of Cooper pairs 
[248, 14], which form a Bose condensate and give rise to an energy gap leading to 
superfluidity, or, in the presence of a gauge symmetry which is spontaneously broken 
by the condensate, to superconductivity. Whereas in conventional superconductors, 
the only attractive interaction is from phonon-mediated interactions, in QCD 
already one-gluon exchange in the colour-antitriplet channel is attractive, so colour- 
superconductivity should be a much less fragile phenomenon with correspondingly 
large energy gaps. This was first studied by Barrois [249] in the late seventies and 
more extensively by Bailin and Love [250] and others [251, 252]. 

Over the last few years there has been a renewed flurry of activity in this 
field following the observation by various groups [253, 254] that the energy gap of 
colour superconductors as well as their critical temperature may be much larger than 
previously expected. Moreover, new symmetry breaking schemes were discovered. 
Whereas it was previously [250] thought that only the lightest quark flavours u and 
d are able to form Cooper pairs, it was argued in reference [253] and corroborated 
in reference [255] that a condensate which locks the breaking of colour and flavour 
symmetry ( "colour-flavour locking" (CFL)) is energetically favoured, at least for high 
baryon chemical potential. In astrophysical applications, where in contrast to heavy- 
ion physics the strange-quark chemical potential is constrained by electric neutrality, 
CFL appears to be even more favoured than two-flavour superconductors (2SC) 
according to the estimates of [256, 257]. Alternatively, electric neutrality may be 
realized with the so-called gapless 2SC phase [258, 259]. 

Close to the phase transition to confined nuclear matter, there is also the 
possibility for rather involved symmetry breaking patterns. A particularly fascinating 
one is the so-called LOFF phase [260, 261, 262], which has a gap function that breaks 
translational invariance, leading to crystalline structures in a colour superconductor 
with possible relevance to neutron star physics [263]. 

Most of the original work has been based on a QCD inspired phenomenological 
Nambu-Jona-Lasinio (NJL) model [264]. Starting with the work of Son [265], the 
phenomenon of colour superconductivity has also been investigated by systematic 
perturbative techniques based on the fundamental Lagrangian of QCD. Son 
established, using renormalization group techniques, that in QCD the parametric form 
of the energy gap is not of the order of b\i cxp(— c/g 2 ) as with point-like interactions, 
but is modified by long-range colour magnetic interactions to the parametrically larger 
order of &i/i<?~ 5 cxp(— Ci/g), with c\ = 37r 2 /\/2. 

Using weak-coupling methods, b\ has been calculated in [266, 267, 268, 269, 270, 
271], and gauge parameter independence of b\ has been verified in [272] in Coulomb- 
like gauges, though covariant gauges present a problem [273, 274, 275, 276]. 

The gauge independence identities of the previous section are not immediately 
applicable, but require a Nambu-Gor'kov ansatz [250, 271, 277] for the inverse 
propagator of the form 



where $ ± are the gap functions, related by <3? (q) = 7o[ < f )+ (q)Y 7o , an d = 




(5.29) 
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Sdet(S^) + ( 5i«°- T _det(5 i ; 1 ) = 6 tot det(S^) = - det(S^)[5X k k + 5X\], (5.31) 



C[S(— g)] T C~ 1 with the charge conjugation matrix C. Flavour and fundamental 
colour indices are suppressed in (5.29). This inverse propagator is the momentum 
space version of the second derivative of the effective action, 

<5§(x)(5*(y) v=^=A?=o,Ag=Ag' 

where * = (ijj, ip c ) T , ^ = (-0,V> C ), and Aq is the expectation value of Aq, which 
is generally nonvanishing in the colour superconducting phase [278]. The doubling 
of fermionic fields in terms of $ and ^ is just a notational convenience here; the 
effective action itself should be viewed as depending only on either (V>,V>) or the set 

*=(V^c) T 

Assuming a spatially homogeneous quark condensate, one can then derive a gauge 
dependence identity for the momentum-space propagator of the form [278] 

d_ 

where the indices i and i comprise colour, flavour, Dirac and Nambu-Gor'kov indices. 
Just like in the case of ordinary spontaneous symmetry breaking, one has to consider 
a total variation [206, 207] of the determinant of the inverse quark propagator, with 
the first term corresponding to the explicit variation of the gauge fixing function, and 
the second term coming from the gauge dependence of Aq. 

Since the determinant is equal to the product of the eigenvalues, equation 
(5.31) implies that the location of the singularities of the quark propagator is gauge 
independent, provided the singularities of 5X k k do not coincide with those of the quark 
propagator. As above, one may argue that SX is 1PI up to a full ghost propagator, 
and up to gluon tadpole insertions, and the singularities of the ghost propagator are 
not correlated to the singularities of the quark propagator. Gauge independence of 
the zeros of the inverse fermion propagator then follows provided that also the 1PI 
parts of SX have no singularities coinciding with the singularities of the propagator. 

At leading order, when the quark self energy can be neglected, this implies that 
the gap function is gauge independent on the quasiparticle mass shell (though at 
higher orders it becomes necessary to consider the complete dispersion relations of the 
quasiparticles). 



6. Hard-thermal-loop effective action and resummation 

As we have seen in several examples now, complete weak-coupling expansions in 
thermal field theory tend to require a reorganization of the standard loop expansion. 
Already the perturbative expansion of the thermal mass (3.8) of a scalar quasi- 
particle in Sect. 3.3 has shown that ordinary perturbation theory fails to determine 
higher-order corrections, but runs into infra-red problems. A resummation of the 
leading-order (HTL) mass is necessary (and sufficient in this case) to reorganize the 
perturbation series [64]. A thermal mass oc gT introduces an additional (soft for 
j < 1) mass scale, and whenever loop calculations receive important contributions 
from this scale, it is clearly mandatory to use propagators dressed by these masses. 

In the general case, it is however equally important to include vertex corrections. 
In gauge theories, this is only natural as Ward identities tie up vertex functions 
with self-energies. But regardless of gauge symmetry considerations, if there are 
contributions to iV-point one-loop vertex functions that are proportional to T 2 like the 
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HTL self-energies, they are as important as bare vertices when the external momentum 
scale is ~ gT. One then has 

r) N C i 

- ' , (6.1) 



T H W TL ~ g N T 2 k 2 - N ~ g N - 2 k 4 - N 



dA N 



for bosonic fields A. 

A one-loop vertex function whose leading contribution (for soft external momenta) 
is proportional to a power of temperature greater than one is called HTL, as this is 
again dominated by a hard loop momentum ~ T. In fact, already in spinor QED 
there are infinitely many HTL vertex functions, namely those involving two external 
fermion lines and an arbitrary number of gauge bosons. In QCD, there are in addition 
HTL's with an arbitrary number of external gluons. These have been first identified in 
references [279, 213, 280] and used to set up a resummation programme for amplitudes 
involving soft external momenta [213, 281]. 

6.1. HTL effective action 

In the case of scalar </> 4 theory, the only HTL is a mass term corresponding to 
a local HTL effective Lagrangian £^ia.r = — i^th^ 2 with m 2 h the leading-order 
term from (3.8). Remarkably, the infinitely many HTL diagrams of gauge theories 
have a comparatively simple and manifestly gauge-invariant integral representation 
[282, 283, 195, 284] 

£HTL _ £HTL _|_ £HTL 

-of dQ v - , v„ rh 2 n f dCl v „„ v n v^ 

= M *J (6 - 2) 

where v = (1, v) is a light-like 4-vector, i.e. with v 2 = 1, and its spatial components 
are averaged over by J dtt v ■ ■ - . Here v is the remnant of the hard plasma constituents' 
momenta ~ Tv^, namely their light-like 4- velocity, and the overall scale T 
has combined with the coupling constant to form the scale of the thermal masses, 
M,m D — gT. 

The covariant derivatives in the denominators of (6.2) are responsible for the 
fact that there are infinitely many HTL's. Because in QED one has D a dj.(A) — > d, 
the only HTL with exclusively photons as external lines is the photon self-energy 
polarization tensor; the other HTL diagrams of QED have two external fermion lines 
and an arbitrary number of photon insertions. 

The gauge boson part of (6.2) has in fact been obtained originally in a form which 
is not obviously gauge invariant, namely 

^ TL =mf,tr|^ + / ^W(vA)} (6.3) 

where the functional W is determined by being a gauge invariant completion of the 
Debye mass term [282] as 

i i °° n 

W(vA) = d (vA)F( — [ivA,*])—vA, Ffzl^V^-. (6.4) 
v ■ o v ■ a ^-^ n + 2 

n— 

This has an interesting interpretation as eikonal of a Chcrn-Simons gauge theory 
[285, 286, 287], where however parity violations and a quantization of the coefficient 
of the action are absent because of the angular average. 
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Explicit representations of HTL vertices for gauge bosons are in fact most 
efficiently obtained from expanding (6.3) rather than the manifestly gauge invariant 
form (6.2). The gauge boson self energy (5.17) is then found to be given by 

, f dn v 

g^og^o - ko 



tl^(k)=m 2 D 



4-7T v ■ k 



(6.5) 



and an n-point vertex function by 

f»i;;;»» ■ ■ • , k n ) = 2g n - 2 m 2 D f ■ • . Vfln \ tr(T a " [!—,[. . . , T a >] ■ ■ •]) 



fc? 1 



" t, ■ t, ■ (*i + k 2 ) ' ' ' v.{k 1 + ... + k n - 2 ) + P™**^ 1 ' (6-6) 
The effective Lagrangian (6.2) can be understood as an effective field theory in 
the sense of Wilson's renormalization group [288]. It arises from integrating out, in 
leading order, the effects of the hard momentum modes of the plasma constituents. 
Because, unlike the case of effective field theories in vacuum field theory, these particles 
are real rather than virtual, they give rise to non-localities corresponding to their free, 
light-like propagation. The HTL vertex functions encoded by (6.2) or (6.3) can thus 
be understood as forward scattering amplitudes for hard (collisionless) particles in soft 
external fields [195], thereby explaining their gauge independence which is less obvious 
on a purely diagrammatic level (there it can be understood through the absence of 
HTL ghost self-energy and vertex functions [192]). 

The action provided by (6.2) or (6.3) is Hermitean only in a Euclidean form. 
After analytic continuation there are cuts corresponding to Landau damping, as 
mentioned above in connection with the two-point functions. In order to obtain the 
analytic continuation relevant for linear response theory, this is better performed on 
the effective equations of motions [289], which are also the primary objects in kinetic 
theory. 



6.2. Kinetic theory approach 

While a kinetic-theory derivation of the HTL propagators generalizing the Abelian 
case [214] to nonabclian Yang-Mills theory has been given already in references [290, 
291, 292, 293], a systematic treatment that includes the HTL vertices has been 
developed only in references [294, 295, 284, 296, 297, 298] (see references [6, 299] 
for two recent comprehensive reviews). In a pure-glue nonabelian gauge theory, the 
effective equations of motion are given by 

[D», F^] a = f a = m 2 D J ^v»W a ( Xl v), (6.7) 
[v D,W(x,v)] a = vE Q (x), (6.8) 

where E is the chromo-electric field strength, and W a (x,x) describes the fluctuations 
of the phase space density of hard gluons. At the expense of introducing W a (x,x) as 
a new soft degree of freedom, this formulation permits a completely local description 
of the physics of soft gauge boson modes [300, 301]. 

The HTL diagrams can be constructed from solving (6.7) in terms of j^[A] and 
formally expanding 

3% = -K A b + \Kv P A b A c +■■■■ ( 6 - 9 ) 
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The kinetic-theory approach allows also to consider strong deviations from 
equilibrium. In references [302, 303, 304] the analogue of the HTL self-energies with 
anisotropic momentum distributions have been considered, which turn out to lead to 
space-like singularities in the gauge boson propagator. The latter are related to plasma 
instabilities which could play an important role in thcrmalization issues [305, 306]. 

6.3. HTL /HDL resummation 

Since for soft momenta <~ gT HTL self-energies and vertices are equally important as 
the tree-level self-energies and vertices, the former may not be treated perturbatively, 
but should rather be combined with the latter to form effective self-energies and 
vertices. This can be done formally by replacing 

£ cl - £ c i + £ HTL - 5 x £ HTL (6.10) 

where 5 is a parameter that is sent to 1 in the end, after the last term has been treated 
as a 'thermal counterterm' by assuming that 5 counts as a one-loop quantity. 

Because £ HTL has been derived under the assumption of soft external momenta, 
this prescription is in fact only to be followed for soft propagators and vertices [213]. 
Propagators and vertices involving hard momenta (if present) do not require this 
resummation, and in fact for obtaining a systematic expansion in the coupling g, one 
has to expand out all HTL insertions on hard internal lines. In practice, a separation 
between hard and soft scales may be implemented by introducing an intermediate scale 
A with jT « A « T, assuming g < 1, for example A oc ^fgT (see e.g. [307, 308]). 

The resulting systematic expansions in g typically involve single powers in g and 
logarithms of g, in contrast to conventional perturbation theory which would involve 
only g 2 as an expansion parameter. This is because increasing the loop order by one 
involves a factor g 2 T which in an ultrarelativistic situation and soft external momenta 
< gT is made dimensionless by a thermal mass m ~ gT, so that the effective expansion 
parameter becomes g 2 T/m <~ g. 

The same resummation scheme arises in the presence of a chemical potential. In 
the effective action, the chemical potential enters only in the mass parameters M 2 
and rh 2 D according to (5.28) and (5.18). For T w but large fif, the HTL have also 
been nicknamed "hard dense loops" (HDL) [309]. At T w 0, resummation of the HDL 
effective action is necessary for soft momenta < gfj,, but because of the absence of 
Bose enhancement this does not give rise to single powers of g in the perturbations 
series, but only to logarithms of g in addition to powers of g 2 , as found long ago in 
the ring resummation scheme of Gell-Mann and Brucckncr for an electron gas at high 
density [310]. HDL resummation for dynamic quantities has been considered first in 
[311] with applications to energy loss rates in astrophysical systems through axion-likc 
particles, and for energy loss of heavy quarks at both large T and [i in [312]. 

HDL resummation is also at the basis of the results of reference [143] quoted 
in Sect. 4.4 on non-Fermi-liquid behaviour of the specific heat at low temperature 
and high chemical potential, equation (4.17). In this case the non-analytic terms in 
g 2 involve a logarithm of g as well as cubic roots of g 2 , which are the result of the 
only weak dynamical screening of near-static magnetic modes whose screening lengths 
involve cubic roots of m 2 D w, see (5.19). 

6.3.1. HTL-screened perturbation theory (HTLPT) In [106, 107, 108, 109] a 
modification of the above scheme has been suggested, where the mass parameters rhp 
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and M within £ HTL are considered as independent of the coupling g and used for a 
variational improvement of the perturbative series. This is a generalization of screened 
perturbation theory [100, 103]. It differs from standard HTL/HDL resummation also 
in that (6.10) is used for both hard and soft momenta, which results in additional UV 
divergences that need to be subtracted at any finite order of the expansion. 

This method has been used to calculate the thermodynamic potential to two loop 
order, where it does improve the apparent convergence of the perturbative results, but 
the result deviates significantly from lattice results even at the highest temperatures 
that are available for the latter. 

A problem of this approach, at least in quantities that are dominated by hard 
excitations like the pressure, seems to be that the HTL effective action is not a 
good approximation at hard momenta. This is also signalled by the fact that HTL 
propagators do not satisfy the relativistic KMS condition at large momenta [17] . While 
this is taken care of eventually by the counterterms in (6.10), at any finite order of 
the expansion there are uncancelled unphysical hard contributions. In [164] it has 
been shown that a (simpler) implementation of a variational perturbation theory in 
dimensional reduction which uses just the Debye mass term (which is the static limit 
of the HTL effective action) avoids this problem and indeed leads to results which 
are closer to the lattice result as well as the higher-order calculations in dimensional 
reduction when improved as discussed in section 4.2. It therefore appears that HTLPT 
needs to be amended such that a different treatment of hard and soft modes is secured. 



6.3.2. HTL resummed thermodynamics through $- derivable approximations While 
HTLPT when applied to the thermodynamical potential does not work satisfactorily 
(though it may be of more use in quantities which depend more dominantly on soft 
rather than hard scales), it turns out that a generalization of the self-consistent 
expression for the entropy (3.18) allows for a resummation of HTL propagators without 
the problems of HTLPT and with remarkably good numerical results when compared 
to lattice data [96, 98, 99]. If 

In gauge theories including fermions, the self-consistent two-loop expression for 
the entropy (3.18) reads 

5=-tr J T^^jS^ [Im log D~\w, k) - ImH(c, k) Re D{u>, k)] 

-2 try 704^^ [lmlogS-V,fc)- ImE(w, fe) ReS(w, fe)] , (6.11) 

and a similarly simple expression can be obtained for the quark number density 

M=-2tvj [Imlogg-V.fc)- ImS(o;, k) Re S(w,k)] . (6.12) 

Here n{u) = (e 0u > - l)" 1 , f{w) = {e^~^ + l)" 1 , and "tr" refers to all discrete labels, 
including spin, colour and flavour when applicable. 

In nonabelian gauge theories, the above expressions have to be augmented by 
Faddeev-Popov ghost contributions which enter like bosonic fields but with opposite 
over-all sign, unless a gauge is used where the ghosts do not propagate such as 
in axial gauges. But because ^-derivable approximations do not generally respect 

% A similar approach but formulated directly in terms of the (modified) thermodynamic potential SI 
has been set up in [313]. 
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gauge invariance, + the self-consistent two-loop approximation will not be gauge- 
fixing independent, ft is in fact not even clear that the corresponding gap equations 
(3.16) have solutions at all or that one can rcnormalizc these (nonperturbative) 
equations, although nonperturbative renormalizability has been proven in the scalar 
case [92, 93, 94, 95]. Concerning gauge fixing dependences it is at least possible to 
show that at a stationary point these enter at twice the order of the truncation [316]. 

In [96, 98, 99] a manifestly gauge invariant approximation to full self-consistency 
has been proposed which maintains equivalence with conventional perturbation theory 
up to order g 3 , which is the maximum (perturbative) accuracy of a two- loop $- 
derivable approximation. For these approximations it will be sufficient to keep only 
the two transverse structure functions of the gluon propagator and to neglect ghosts. 

For soft momenta, the appropriate leading order propagator is the HTL one, and 
indeed there is no HTL ghost self-energy. 

For hard momenta, one can identify the contributions to (6.11) below order g 4 as 
those linear in the self-energies, 

— -j -J- (ReS+Im— - - Rc£_Im — -} (6.13) 
(2tt) 4 dT\ + Lu-k Lu + ki v ' 

considering now a gauge theory with N g gluons and Nf fermion flavours. Because the 
imaginary parts of the free propagators restrict their contribution to the light-cone, 
only the light-cone projections of the self-energies enter. At order g 2 this is exactly 
given by the HTL results, without having to assume soft w, k [317, 318] 

Renf \lo 2 = k 2 ) = n t (w 2 = k 2 ) = \m\ = m^, (6.14) 

2k Rc54 2) (w = ±fc) = 2fcS±(w = ±k) = 2M 2 = M^, (6.15) 

and without contributions from the other components of and the Faddeev-Popov 
self-energy. 

There is no contribution oc g 2 from soft momenta in (6.11) and (6.12) so that 
one is left with remarkably simple general formulae for the leading-order interaction 
contributions to the thermodynamic potentials expressed through the asymptotic 
thermal masses of the bosonic and fermionic quasiparticles: 

^ {2) =-^J2^ M2 ooF- (6-16) 

F 

Here the sums run over all the bosonic (B) and fermionic (F) degrees of freedom (e.g. 
4 for each Dirac fermion), which are allowed to have different asymptotic masses and, 
in the case of fermions, different chemical potentials. 

The result (6.16) also makes it clear that relative-order-^ corrections to m 2 ^, 
Mx> ~ 9 2 T 2 will contribute to the order-g 3 terms in S and N. 

The HTL approximation to S and TV thus includes correctly the leading-order 
interaction term oc g 2 and only part of the order-g 3 terms. Using the peculiar 
sum rule (A. 13) one can in fact show that in the case of pure-glue QCD the HTL 
entropy contains exactly 1/4 of the plasmon term ~ g 3 . When treated strictly 
perturbatively, even 1/4 of the plasmon term spoils the apparent convergence. 

+ For this, one would have to treat vertices on an equal footing with self-energies, which is in principle 
possible using the formalism developed in references [314, 139, 315]. 
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However, in the nonperturbative expression (6.11) the otherwise large g 3 correction 
is rendered harmless and leads to a small correction such that the rough agreement 
of the perturbative order-g 2 result with lattice results for T > 3T C is retained and 
improved. 

The plasmon term ~ g 3 becomes complete only upon inclusion of the next- 
to-leading correction to the asymptotic thermal masses moo and Moo- These are 
determined in standard HTL perturbation theory through 

5m£,(fc) - Re SU t (lu = k) 

-l j*^ ( 6 - 17 ) 

= Re( 

where thick dashed and wiggly lines with a blob represent HTL propagators for 
longitudinal and transverse polarizations, respectively. Similarly, 

^SMl(k) = 5Z+(gj = k) = Re( , lTX + . )\ u=k . (6.18) 

The explicit proof that these contributions indeed restore the correct plasmon term is 
given in reference [99]. 

These corrections to the asymptotic thermal masses are, in contrast to the latter, 
nontrivial functions of the momentum, which can be evaluated only numerically. 
However, as far as the generation of the plasmon term is concerned, these functions 
contribute in the averaged form 

= Jdkkn' (k)ReSU T (^k) 
°° Jdkkn' BE (k) 
(cf. (6.13)) and similarly 

- 2 _ Jdkkn' FD (k)Re2kSJ: + (u J = k) 

6M °°~ Jdkkn' FB (k) ■ (b - M) 

These averaged asymptotic thermal masses turn out to be given by the remarkably 
simple expressions [99] 

5m 2 x = -±-g 2 NTm D , SM^ = -±-g*C f Tm D , (6.21) 

where C/ = N g /(2N). Since the integrals in (6.19) and (6.20) are dominated by hard 
momenta, these thermal mass corrections only pertain to hard excitations. 

Pending a full evaluation of the NLO corrections to Re^n and Re5£, 
references [96, 98, 99] have proposed to define a next-to-leading approximation through 
(for gluons) 



Snla — Shtl 



soft + s HTL , m ^ 



(6.22) 

hard 



where m 2 ^ includes (6.21) and a separation scale \fc^2-KTmD is introduced to make 
the distinction between hard and soft domains. 

For numerical evaluations, a crucial issue here is the definition of the corrected 
asymptotic mass moo. For the range of coupling constants of interest (g > 1), the 
correction l^m^l is greater than the LO value m 2 ^, leading to tachyonic masses if 
included in a strictly perturbative manner. 

However, this problem is not at all specific to QCD. In section 3.5 we have seen 
that the perturbative result (3.8) for the scalar screening mass to order g 3 also turns 
tachyonic for g > 1. The self-consistent one-loop gap equation (3.6), on the other 
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Figure 11. Comparison of the lattice data for the entropy of pure-glue SU(3) 
gauge theory of reference [154] (gray band) with the range of Shtl (solid lines) 
and Smla (dash-dotted lines) for p, = nT . . . 4nT and ca = 1/2 ... 2. 

hand, is a monotonic function in g, and is well approximated by a quadratic equation 
(3.13) which keeps just to first two terms in in m/T expansion. As a model for the 
rather intractable (nonlocal) gap equations of QCD in a ^-derivable approximation, 
reference [99] has therefore proposed to include the corrections (6.21) in analogy to 
(3.13) in order to enforce monotonicity of TTIqq clS function of the coupling. (In order 
to be consistent with the large flavour number limit, this should in fact be done only 
for the bosonic part; for Nf ^ 0, the fermionic gap equation has to remain linear in 
M%o, which however poses no problem at small Nf [319].) 

In figure 11, the numerical results for the pure-glue HTL entropy and the NLA 
one (using a quadratic gap equation for moo) are given as a function of T/T c with 
T c chosen as T c = 1.14A^jg. The full lines show the range of results for Shtl when 
the renormalization scale p, is varied from nT to 47rT; the dash-dotted lines mark the 
corresponding results for Snla with the additional variation of ca from 1/2 to 2. The 
dark-gray band are lattice data from reference [154]. Evidently, there is very good 
agreement for T > 2.5T C . 

This approach can be generalized [98, 99, 171] also to nonzero chemical potentials 
Hf, for which lattice results are available in the form of quark number susceptibilities 
[320, 321, 322, 323, 324] and are beginning to become available for finite small 
chemical potential [173, 174, 175, 176, 177, 178, 179, 180]. Simpler quasiparticle 
models [160, 161] have already been used to extrapolate lattice data to finite chemical 
potential [162] and seem to work well when compared with the recent lattice results 
for nonzero chemical potential [325]. The HTL approach offers a possible refinement, 
which has been worked out in [326, 327]. 

7. Next-to-leading order corrections to the quasi-particle spectrum 

In the simple scalar model of Sect. 3.3 and in the more complicated example of 
asymptotic (averaged) gluon and quark masses, eqs. (6.19) and (6.20), we have seen 
that HTL/HDL perturbation theory leads to next-to-leading order corrections of quasi- 
particle dispersion laws which are typically suppressed by a single power of g rather 
than g 2 . 
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In the following, we shall review what is at present known about next-to-leading 
order corrections to the HTL quasi-particles in gauge theories. It will turn out that 
some corrections, namely screening lengths for frequencies below the plasma frequency 
and damping rates for moving excitations are even more enhanced by infrared effects, 
to an extent that they cannot be determined beyond the leading logarithmic term 
without nonperturbative input. 

7.1. Long-wavelength plasmon damping 

Historically, the first full-fledged application of HTL resummation techniques was the 
calculation of the damping constant of gluonic plasmons in the long- wavelength limit. 
In fact, the development of the HTL resummation methods was stimulated by the 
failure of conventional thermal perturbation theory to determine this quantity. 

In bare perturbation theory, the one-loop long-wavelength plasmon damping 
constant had been found to be gauge-parameter dependent and negative definite 
in covariant gauges [216]. Since even the gauge-independent frameworks of the 
Vilkovisky-DeWitt effective action [328, 329] and the gauge-independent pinch 
technique [330] led to negative one-loop damping constants, this was sometimes 
interpreted as a signal of an instability of the quark-gluon plasma [331, 332]. Other 
authors instead argued in favour of ghost-free "physical" gauges such as temporal 
and Coulomb gauge, where the result turned out to be positive and seemingly gauge- 
independent in this class of gauges [196, 197]. 

It was pointed out in particular by Pisarski [333] that all these results at bare one- 
loop order were incomplete, as also implied by the arguments for gauge independence 
of the singularities of the gluon propagator of reference [191]. The appropriate 
resummation scheme was finally developed by Braaten and Pisarski in 1990 [281, 213] 
who first obtained the complete leading term in the plasmon damping constant by 
evaluating the diagrams displayed in figure 12 [334] with the result 

7 (k = 0) = -?-Im 6U A (fc = u pU k = 0) 

ZLUpl 

= +6.635...^-^ = 0.264^5 w p i. (7.1) 

247T 

For g <C 1 this implies the existence of weakly damped gluonic plasmons. In QCD 
(N — 3), where for all temperatures of interest g > 1, the existence of plasmons as 
identifiable quasi-particles requires that g is significantly less than about 2.2, so that 
the situation is somewhat marginal. 

The corresponding quantity for fermionic quasi-particles has been calculated in 
[335, 336] with a comparable result: weakly damped long-wavelength fermionic quasi- 
particles in 2- or 3-flavour QCD require that g is significantly less than about 2.7. 

These results have been obtained in Coulomb gauge with formal verification 
of their gauge independence. Actual calculations, however, later revealed that in 
covariant gauges HTL-resummed perturbation theory still leads to explicit gauge 
dependent contributions to the damping of fermionic [210] as well as gluonic [211] 
quasi-particles. But, as was pointed out subsequently in [212], these apparent gauge 
dependences are avoided if the quasi-particle mass-shell is approached with a general 
infrared cut-off (such as finite volume) in place, and this cut-off lifted only in the 
end. This procedure defines gauge-independent dispersion laws; the gauge dependent 
parts are found to pertain to the residue, which at finite temperature happens to be 
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Figure 12. One-loop diagrams in HTL-rcsummcd perturbation theory 
contributing to the gauge boson self-energy in pure-glue QCD. HTL-rcsummcd 
quantities are marked with a blob. 



linearly infrared singular in covariant gauges, rather than only logarithmically as at 
zero temperature, due to Bose enhancement. 

For kinematical reasons, the result (7.1), which has been derived for branch A 
of the gluon propagator, should equally hold for branch B, since with k — > one 
can no longer distinguish between spatially transverse and longitudinal polarizations. 
However, as will be discussed further below, the limit k — > involves infrared problems, 
and there are even explicit calculations [337, 338] that claim to find obstructions to 
this equality, which are however contradicted by the recent work of [339]. 

In the simpler </> 4 theory the long-wavelength plasmon damping constant has been 
calculated in [64, 65] and it has been shown in [340, 341, 342, 343] that this quantity 
can also be extracted from classical field theory after perturbativcly matching to the 
HTL mass. 

7.2. NLO correction to the plasma frequency 

In [344] , Schulz has calculated also the real part of the NLO contribution to the gluon 
polarization tensor in the limit of k — > which determines the NLO correction to the 
gluonic plasma frequency. 

The original power-counting arguments of [213] suggested that besides one-loop 
diagrams with HTL-resummed propagators and vertices given in figure 12, there could 
be also contributions from two-loop diagrams to relative order g. The explicit (and 
lengthy) calculation of [344] showed that those contribute only at order g 2 ln(l/g) 
rather than <?, and the NLO plasma frequency in a pure-glue plasma was obtained as 



In this particular result, HTL-resummed perturbation theory turns out to give a 
moderate correction to the leading-order HTL value even for g <~ 1. 

7.3. NLO correction to the Debye mass 

Poles of the gauge boson propagator at lu < u p \. and k 2 = — k 2 < describe 
exponential screening of fields with frequencies below the plasma frequency and, for 
lu — > 0, of static fields. The value of n is the (frequency-dependent) screening mass. 

In the static case, branch A of the propagator describes the screening of (chromo-) 
magnetostatic fields. While there is a finite screening length as long as lu > 0, the 
A-branch of the HTL propagator becomes unscreened in the static limit. In QED, a 
"magnetic mass" is forbidden by gauge invariance [215, 222], but some sort of entirely 
non-pcrturbative magnetic mass is expected in non-Abclian gauge theories in view 
of severe infrared problems caused by the self-interactions of magnetostatic gluons 



oJpi. = oJpi. 1 - 0.09VNg . 



(7.2) 



[149, 150, 151]. 
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Branch B, on the other hand, contains the information about screening of 
(chromo) electric fields as generated by static charges (Dcbyc screening). The Debye 
mass given by the leading-order HTL propagator is m D = V3lj p \. . 

The evaluation of corrections in thermal perturbation theory again requires 
resummation. Even though the nonlocalities of the HTL effective action do not play a 
role in the static limit, historically these corrections have been determined for QCD, 
as well as for (ultra-relativistic) QED, only after the more complicated cases of gluonic 
plasmon damping and NLO plasma frequency discussed above had been mastered. 

Originally, the Dcbyc mass (squared) has been defined as the infrared limit 
n o(w = 0, k — > 0), which indeed is correct at the HTL level, cf. (5.17). 

In QED, this definition has the advantage of being directly related to the second 
derivative of the thermodynamic pressure with respect to the chemical potential 
/j, (electric susceptibility), so that the higher-order terms known from the latter 
determine those of U^ b (uj = 0, k -> 0) through [215, 8] 



dfi 2 



This result is gauge independent because in QED all of is. 

In the case of QCD, there is no such relation. Moreover, 5m 2 D /m 2 D ~ g rather 
than g 3 because of gluonic self-interactions and Bose enhancement. The calculation of 
this quantity should be much easier than the dynamic ones considered above, because 
in the static limit the HTL effective action collapses to just the local, bilinear HTL 
Debye mass term, 

£ HTL ^tic _^2 ntrA l {1A) 

This is also gauge invariant, because A$ behaves like an adjoint scalar under 
time-independent gauge transformations. Resummed perturbation theory for static 
quantities thus reduces to a resummation of the HTL Debye mass in the electrostatic 
propagator [310, 8, 114]. 

However, in QCD this simple ("ring") resummation leads to the gauge dependent 
result [345] 




n«K,(o.())Mi) = i -'-'-tJ^ftV (7 ' 5) 

where a is the gauge parameter of general covariant gauge (which coincides with 
general Coulomb gauge in the static limit). 

This result was initially interpreted as meaning either that the non-Abelian 
Debye mass could not be obtained in resummed perturbation theory [346] or that 
one should use a physical gauge instead [196, 8]. In particular, temporal axial gauge 
was put forward, because in this gauge there is, like in QED, a linear relationship 
between electric field strength correlators and the gauge propagator. However, because 
static ring resummation clashes with temporal gauge, inconclusive and contradicting 
results were obtained by different authors [347, 348, 196]. A consistent calculation 
in fact requires vertex resummations [349, 350], but this does not resolve the gauge 
dependence issue because the nonabelian field strength correlator is gauge variant 
[351]. 

On the other hand, in view of the gauge dependence identities discussed in section 
5.1.1, the gauge dependence of (7.5) is no longer surprising. Gauge independence can 
only be expected "on-shell", which here means u = but k 2 — > —m 2 D . 
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Indeed, the exponential fall-off of the electrostatic propagator is determined by 
the position of the singularity of Ag(0, k), and not simply by its infrared limit. This 
implies in particular that one should use a different definition of the Debye mass 
already in QED, despite the gauge independence of (7.3), namely [352] 

m 2 D = n 0O (0,k) (7.6) 

For QED (with massless electrons), the Debye mass is thus not given by (7.3) but 
rather as [352] 

m 2 D = n 00 (o,fc^o) + [n 00 (o,A;)| fe2= _ mi) -n 00 (o,fc^o)] 

e 2 T 2 / 3e 2 V3e 3 e 2 n Jl 4. \ 

where jl is the renormalization scale of the momentum subtraction scheme i.e. 
n M „(fc 2 = — /1 2 )|t=o = 0. In the also widely used modified minimal subtraction (MS) 
scheme the last coefficient — | in (7.7) has to be replaced by — i. (The slightly different 
numbers in the terms oc e 4 T 2 quoted in [222, 10] pertain to the minimal subtraction 
(MS) scheme*.) 

Since de/d\nfl = e 3 /(127r 2 ) + 0(e 5 ), (7.7) is a renormalization-group invariant 
result for the Debye mass in hot QED, which (7.3) obviously is not. Only the 
susceptibility \ = d 2 P/d/j, 2 is renormalization-group invariant, but not e 2 %- Similarly, 
one should distinguish between electric susceptibility and electric Debye mass also in 
the context of QCD. 

In QCD, where gauge independence is not automatic, the dependence on the 
gauge fixing parameter a is another indication that (7.5) is the wrong definition. 
For (7.6) the full momentum dependence of the correction SU no (ko — 0, k) to IIoo is 
needed. This is given by [352] 



<m 00 (fco = 0, k) = gm D NJ m 6 + — J 



6 f d 3 ~ 2e p ( 1 1 
~Nf J (2tt) 3 - 2£ 1 p 2 + m 2 D + p 2 



D 

4m 2 3 -(k 2 +m 2 3 )[3 + 2pk/p 2 ] 2 2 p 2 + 2pk 

p 2 [(p + k) 2 + m 2 D ] +a[k +mcj p4[(p + k) 2 +™ 2 D ]. 

In accordance with the gauge dependence identities, the last term shows that 
gauge independence holds algebraically for k 2 = —m 2 D . On the other hand, on this 
"screening mass shell", where the denominator term [(p + k) 2 + mf,] — ► [p 2 + 2pk], 
one encounters IR-singularities. In the a-dependent term, they are such that they 
produce a factor 1/ [k 2 + m 2 D ] so that the gauge dependences no longer disappear even 
on-shell. This is, however, the very same problem that had to be solved in the above 
case of the plasmon damping in covariant gauges. Introducing a temporary infrared 
cut-off (e.g., finite volume), does not modify the factor [k 2 -!-^. 2 }] in the numerator but 
removes the infrared divergences that would otherwise cancel it. Gauge independence 
thus holds for all values of this cut-off, which can be sent to zero in the end. The gauge 
dependences are thereby identified as belonging to the (infrared divergent) residue. 

The third term in the curly brackets, however, remains logarithmically singular 
on-shell when the infrared cut-off is removed. In contrast to the a-dependent 
term, closer inspection reveals that these singularities are coming from the massless 
magnetostatic modes and not from unphysical massless gauge modes. 

* Reference [10] erroneously refers to the MS scheme when quoting the MS result of [222] . 
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At HTL level, there is no (chromo) magnetostatic screening, but, as we have 
mentioned, one expects some sort of such screening to be generated non-perturbatively 
in the static sector of hot QCD at the scale g 2 T <~ gm D [149, 150, 151]. 

While this singularity prevents evaluating <5n o in full, the fact that this 
singularity is only logarithmic allows one to extract the leading term of (7.8) under 
the assumption of an effective cut-off at p ~ g 2 T as [352] 



The (9((7)-contribution, however, is sensitive to the physics of the magnetostatic 
sector at scale g 2 T, and is completely non-pcrturbative in that all loop orders > 2 are 
expected to contribute with equal importance as we shall discuss in section 7.4. 

Because of the undetermined 0(g)-term in (7.9), one-loop resummed perturbation 
theory only says that for sufficiently small g, where 0(g\n(l/g)) >• 0(g), there is a 
positive correction to the Debye mass of lowest-order perturbation theory following 
from the pole definition (7.6), and that it is gauge independent. 

On the lattice, the static gluon propagator of pure SU(2) gauge theory at high 
temperature has been studied in various gauges [353, 223] with the result that the 
electrostatic propagator is exponentially screened with a screening mass that indeed 
appears to be gauge independent and which is about 60% larger than the leading- 
order Debye mass for temperatures T/T c up to about 10 4 . Similar results have been 
obtained recently also for the case of SU(3) [354, 355]. 

In [351], an estimate of the 0(g) contribution to (7.9) has been made using the 
crude approximation of a simple massive propagator for the magnetostatic one, which 
leads to 



On the lattice one finds strong gauge dependences of the magnetostatic screening 
function, but the data are consistent with an over-all exponential behaviour 
corresponding to m m w 0.5g 2 T in all gauges [353, 356]. Using this number in a self- 
consistent evaluation of (7.10) gives an estimate for m D which is about 20% larger 
than the leading-order value for T/T c = 10 . . . 10 4 . 

This shows that there are strong non-perturbative contributions to the Debye 
screening mass mo even at very high temperatures. Assuming that these are 
predominantly of order g 2 T, one-loop resummed perturbation theory (which is as far as 
one can get) is able to account for about 1 /3 of this inherently non-perturbative physics 
already, if one introduces a simple, purely phenomenological magnetic screening mass. 

7.3.1. Non-perturbative definitions of the Debye mass A different approach to 
studying Debye screening non-perturbatively without the complication of gauge fixing 
is to consider spatial correlation functions of appropriate gauge-invariant operators 
such as those of the Polyakov loop 



The correlation of two such operators is related to the free energy of a quark-antiquark 
pair [357]. In lowest order perturbation theory this is given by the square of a Yukawa 




(7.9) 




(7.10) 




(7.11) 
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potential with screening mass rh D [346]; at one- loop order one can in fact identify 
contributions of the form (7.10) if one assumes magnetic screening [358, 351], but 
there is the problem that through higher loop orders the large-distance behaviour 
becomes dominated by the magnetostatic modes and their lightest bound states [359]. 

In [360], Arnold and Yaffe have proposed to use Euclidean time reflection 
symmetry to distinguish electric and magnetic contributions to screening, and have 
given a prescription to compute the sublogarithmic contribution of order g 2 T to mo 
nonperturbativcly. This has been carried out in 3-d lattice simulations for SU(2) 
[361, 362] as well as for SU(3) [363]. The Debye mass thus defined shows even larger 
deviations from the lowest-order perturbative results than that from gauge- fixed lattice 
propagators. E.g., in SU(2) at T = 10 4 T C this deviation turns out to be over 100%, 
while in SU(3) the dominance of g 2 T contributions is even more pronounced. 

Clearly, (resummed) perturbation theory is of no use here for any temperature of 
practical interest. However, it should be noted that the magnitude of the contributions 
from the completely nonperturbative magnetostatic sector depends strongly on the 
quantity considered. It is significantly smaller in the definition of the Debye mass 
through the exponential decay of gauge- fixed gluon propagators, which leads to smaller 
screening masses on the lattice [353, 223, 354, 355] and which also seem to provide a 
useful nonperturbative description of Debye screening as they turn out to be gauge- 
independent in accordance with the arguments of section 5.1.1. 

7.4- Magnetostatic screening 

The next-to-leading order correction to the other structure function, ILi(0, k), can be 
derived in a calculation analogous to the one leading to (7.8), and is found to be [352] 



Apparently, the magnetic permeability defined by = 1 — \\\u/k 2 is a gauge- 
dependent quantity beyond leading order. The gauge-dependent terms vanish only at 
the location of the pole of the transverse gluon propagator, which is at k = 0. There 
the correction term vanishes completely which means that there is no magnetic mass 
squared of the order g 3 T 2 . The magnetic mass must therefore be <C g 3 ^ 2 T at weak 
coupling. 

For small k <C Ad, (7.12) has a linear behaviour with gauge dependent, 
but positive definite coefficient. The transverse propagator therefore has the form 
l/(k 2 — ck) with c ~ g 2 T. This corresponds to a pole at space-like momentum with 
k <~ g 2 T [364]. However, in this regime the inherently nonperturbative contributions 
~ g A T 2 to n^j, become relevant and are expected to remove this pathology by the 
generation of a magnetic mass m m ~ g 2 T. This is completely analogous to the 
infrared problem of the perturbative expansion of the pressure identified in [150, 151], 
but while the latter sets in at 4-loop order, in the transverse gauge-boson self-energy 
with external momenta < g 2 T this problem starts already at 2-loop order. 

In a HTL-resummed 2-loop calculation in general covariant gauge, reference [365] 
has verified that such a magnetic mass receives contributions exclusively from the 
ultrasoft momentum regime k ~ g 2 T where the relevant effective theory is three- 
dimensional Yang-Mills theory, which is confining and inherently nonperturbative. 




(7.12) 
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There exists in fact an elegant attempt towards an analytic nonperturbative study 
of this theory through its Schrodinger equation using results from two-dimensional 
conformal field theory [366, 367, 368]. This leads to the estimate m m — g 2 NT/(2n), 
which is reasonably close to the results for magnetostatic propagators obtained in the 
lattice calculations of references [353, 223] quoted above. 

The other approaches that have been tried to obtain an analytic estimate for the 
magnetic mass mostly involve a reorganization of perturbation theory by assuming 
a more or less complicated mass term for the magnetostatic gluons, setting up a 
self-consistent gap equation and solving it. However, the various possibilities lead to 
rather contradictory results [369, 370, 371, 372, 373, 374]. In fact, lattice calculations 
indicate that the phenomenon of magnetic screening is not well described by a simple 
pole mass [223]. 

In the Abelian case, one may expect that magnetostatic fields are completely 
unscreened and it can indeed be proved rigorously to all orders of perturbation theory 
[215, 222] that 11^(0, A; — > 0) = 0(k 2 ). In massless scalar electrodynamics, an 
unresummed one-loop calculation actually gives Hu(0,k — > 0) = |e 2 fcT suggesting 
that the magnetostatic propagator has the form l/[fc 2 + e 2 kT/8], which would imply 
power-law screening. However, a resummation of the thermal mass of the scalar 
particles removes this unphysical behaviour [375, 222]. 



7.5. Dynamical screening and damping at high temperature 

A logarithmic sensitivity to the nonperturbative physics of the (chromo-) magneto- 
static sector has in fact been encountered early on also in the calculation of the 
damping rate for a heavy fermion [333], and more generally of hard particles 
[376, 377, 378, 379]. It also turns out to occur for soft quasi-particles as soon as 
they have nonvanishing (group) velocity [380, 381]. 

Because this logarithmic sensitivity arises only if one internal line of (resummed) 
one-loop diagrams is static, the coefficient of the resulting g ln(l/<jr)-term is almost 
as easy to obtain as in the case of the Debye mass, even though the external line is 
non-static and soft, therefore requiring HTL-rcsummed vertices (see figure 12). 

The infrared singularity arises (again) from the dressed one-loop diagram with 
two propagators, one of which is magnetostatic and thus massless in the HTL 
approximation, and the other of the same type as the external one, so only the first 
diagram in figure 12 is relevant. The dressed 3-vertices in it are needed only in the 
limit of one leg being magnetostatic and having zero momentum. Because of the 
gauge invariancc of HTL's, these are determined by the HTL self-energies through a 
differential Ward identity, e.g. 

d_ 

dk~e 

for the 3-gluon vertex (colour indices omitted). 
Comparatively simple algebra gives [381] 

aili(fc) -g 2 N4k 2 [l + d k 2ll I (k)} 2 S I (k), I = A,B (7.14) 

where 

f d 3 p 1 -1 

<Sj(*0 : = T J J^j3p2 (fc-p)2_ n/ (fc-p) fe2= nj (fe),p0 = o (7 ' 15) 

and the logarithmic (mass-shell) singularity arises because (k — p) 2 — Hi(k — p) — > 
-p 2 + 2pk -IL I (k-p) + Ilj(fc) - |p| as k 2 Ilj(fc). 



t^ e (k;-k;0) = -— % v {k) (7.13) 
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The IR-singular part of iSj-(fc) is given by 

d 3 p 1 



5/(fc) T J (27r) 3 p 2 p 2 -2pk + n / (fc-p)-n / (fc)-ie 

_rp + 4,n,(t)]-r* i ta £±^ (7.16, 

J x p 2|k| p - 2|k| - le 

where in the last line we have inserted an IR cutoff A <§C gT for the p-integral in order 
to isolate the singular behaviour. 

One finds that (7.16) has a singular imaginary part for propagating modes, and 
a singular real part in screening situations 

The case k = 0, on the other hand, is IR-safe, because (7.14) is proportional to k 2 , 
while 



T 



k=o V A 2 



Tiki 



5/(fc)— »^[l + ft ca n 7 (fc)]- 1 i n + 0(-£-) fork^O. (7.18) 



This shows that there is a common origin for the infrared sensitivity of screening 
and damping of HTL quasi-particlcs. Provided that the scale where the logarithmic 
divergences are cut off is the magnetic scale g 2 T, the coefficients of the leading 
g In (1/g) -terms are determined and in fact lead to beautifully simple results: For 
the damping of moving quasi-particles one obtains [380, 381] 

/„ ,x 9 2nt Iklfl + d^IIjffc)] , 1 g 2 NT ,„, N1 1 

7i k ' k ^ I ., 111-= y ——v I k In- 7.19 

4tt w(|k|)[l - 5 w 2n/(fc)J g 4vr 3 

where t»j(|k|) is the group velocity of mode I (which vanishes at k = 0). The IR- 
sensitive NLO correction to screening takes its simplest form when formulated as 
[381] 

S K j(uj) = L^I. Kl (u) (in- + 0(1)] (7.20) 
2?r V 5 / 

where K/(w) is the inverse screening length of mode I at frequency w < w p i. (which 

in the static limit approaches the Debye mass and perturbatively vanishing magnetic 

mass, resp., while approaching zero for both modes asw^ <*V.). 

A completely analogous calculation for the fermionic modes (for which there are 

no screening masses) gives 

7±(|k|) = ^g^>±(|k|)| (lnl+0(l)) (7.21) 

for |k| > 0. The group velocity v± equals ±^ in the limit k — > 0, and increases 
monotonically towards +1 for larger momenta (with a zero for the (— )-branch at 
|kdip|/M w 0.41). For strictly |k| = 0, the IR sensitivity in fact disappears because 
(7.21) is no longer valid for |k| <C A, but one has 7±(|k|)| s i ng . oc 3 2 T|k|/A instead. 
Thus 7±(0) is calculable at order g 2 T in HTL-resummed perturbation theory, and 
has been calculated in [335, 336]. For nonvanishing |k| ~ A ~ g 2 T HTL-resummed 
perturbation theory breaks down (the IR sensitivity of 7± for small nonvanishing k 
has recently been displayed also in [382]). For |k| > gT, the damping is calculable in 
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Figure 13. Ncxt-to-leading order contributions to the self-energy of a hard 
fcrmion. The dashed and wiggly lines with a blob refer to the HTL/HDL 
resummed longitudinal and transverse gauge boson propagators, respectively. 

the leading-log approximation, and even completely at |k| = |kdi p | in the case of 7_, 
though this has not been done yet. 

The fermionic result (7.21) applies in fact equally to QED, for which one just 
needs to replace g 2 CF — > e 2 . This is somewhat disturbing as QED docs not allow a 
non-zero magnetic mass as IR cutoff, and it has been conjectured that the damping 
7 - g 2 T or e 2 T itself might act as an effective IR cutoff [376, 377, 380, 383], which 
however led to further difficulties [384, 385]. The solution for QED was finally found 
in references [386, 387, 388, 389] where it was shown that there the fermionic modes 
undergo over-exponential damping in the form e -7 * — > exp(— j^Tt ln(wpi.i)) (for 
v — > 1), so finite time is the actual IR cut-off (see also [390, 391, 392]). The fermion 
propagator has in fact no simple quasi-particle pole in momentum space or any other 
singularity near the light-cone [190] , but nevertheless a sharply peaked spectral density. 

In non-Abclian theories, on the other hand, one does expect static (chromo-)- 
magnetic fields to have finite range, and lattice results do confirm this, so the above 
leading-log results for an (exponential) damping constant are expected to be applicable 
after all, at least for sufficiently weak coupling. 

As will be discussed further in section 8, the damping rate of hard gluons 
7 = 7A(|k| <~ T) ~ -^g 2 NT\n(l/ g) defines an important scale in the dynamics of 
nonabelian fields that is parametrically larger than the magnetic mass scale g 2 T. It 
sets the time scale for colour relaxation and determines, to leading order, the colour 
conductivity through [393, 394, 395] 

a c = £ft/7- (7.22) 

7. 6. Damping of high-momentum fermions in a degenerate plasma 

At zero temperature and nonvanishing chemical potential, the damping rate of 
fermionic excitations is calculable in perturbation theory, for in the absence of Bosc 
enhancement the dynamical screening of the quasi-static transverse modes in the HDL 
gauge boson propagator is sufficient to remove all infrared divergences [396, 397, 398]. 
The leading-order damping rate of a high energy /momentum mode with E,p^> g[i is 
given by the diagrams in figure 13, where only the gauge boson propagator needs to 
be resummed. Explicitly, one has 

7+ (£) = ^ J -J^ [e(qo) -e( P -E + qo )] 

x [Pi(qo,q)+v 2 (l-cos 2 0pt(qo,q)] qo=qvcosi (7.23) 

where v = p/E is the velocity of the fermion and p t ,£ are defined in (A.l) and (A. 2). 
(In the case of QED g 2 Cp reduces simply to e 2 .) For massless fermions v = 1 and 
Pt.i are the HDL expressions given in appendix A.l, but (7.23) remains applicable 
to massive fermions and correspondingly modified spectral functions for the gauge 
bosons [309, 398]. 
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7+(£) = y -7^r\E - Ml + - ^ + 0(1^ - Ml 3 ) ( 7 - 24 ) 



Since dynamical screening vanishes in the zero-frequency limit according to (5.19), 
there is a qualitative deviation from the usual behaviour of the damping of fcrmionic 
excitations in the vicinity of the Fermi surface, which for short-range interactions 
vanishing like (E — fi) 2 as E — > \i due to phase-space restrictions in the fermion- 
fcrmion scattering [399]. The magnetic interactions instead give rise to a \E — fi\ 
behaviour characteristic of non-Fermi liquids [183, 396, 397, 398] which reads 

"—?-\E— I I g2< -" F 
24tt 1 Ml 6Av 2 F m D ' 

with m 2 D = vpm 2 D and m 2 D given by the T = limit of (5.18). Here the first term 
is from the dynamically screened quasi-static magnetic interactions and the second 
term from Debye screened electrostatic ones. In a nonrelativistic situation (v <C 1) 
the latter is the dominant one (except very close to the Fermi surface). 
Far away from the Fermi surface, the damping approaches [397] 

-y + (E) ~ 0.019 g 2 C F m D , £-/i>/z (7.25) 

which is also the damping relevant for hard anti-fermions [398] (which are auto- 
matically far from their Fermi surface). 

7. 7. Non-Fermi-liquid contributions to the real part of the fermion self-energy 

The real part of the fermion self-energy at zero temperature can equally be calculated 
in HDL perturbation theory, but the most interesting aspect of it, a nonanalytic 
behaviour in the vicinity of the Fermi surface, can be inferred from the behaviour 
of its imaginary part through a Kramers-Kronig dispersion relation, which implies 
[183, 270, 398, 400] 

Rc£ + (i?,p) ~ RcS+( M ,p) - ^(E - n) In + 0(\E - M |). (7.26) 

This quantity is gauge independent on the mass shell of the hard particle, to which 
(7.26) presents the leading correction. The energy-independent part is governed by 
the asymptotic fermionic mass, 

ReS + (/x,p) = Mi/(2p), M 2 oo = 2M 2 = ^Lg 2 C f ^ 2 . (7.27) 

This corresponds to a correction to the Fermi momentum defined by u + (p F ) = fi, 
which for effectively massless fermions reads 

8ir 2 

The logarithmic term in (7.26) leads to a correction to the group velocity of the form 
dE , g 2 C t - m D 

v ^^- 1 -^Je-JA + — (7 ' 29) 

which dominates over the contribution from the asymptotic fermion mass for \E — 
/Lt| <C tud, and eventually spoils (HDL) perturbation theory when \E — ji\/mD ^ 
exp(— C/g 2 ). Evidently, the quasistatic magnetic interactions lead to significant 
changes in the vicinity of the Fermi surface. Such non-Fermi-liquid corrections can 
even be the dominant effects in certain quantities, such as the entropy or specific 
heat at low temperature T <C bid ~ gfj,, as discussed in section 4.4. In a colour 
superconductor, where the quasiparticle and quasiholes at the Fermi surface develop a 



P^ = l-y + •■•• (7-28) 
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gap of order bi[ig~ 5 exp(— C\/g) as mentioned in section 5.3, the non-Fermi-liquid 
corrections to the real part of the fermion self-energy remain perturbative and 
contribute at the level of the constant 61, resulting in a significant reduction of the 
magnitude of the gap [269, 271] compared to previous results which did not include 
non- Fermi-liquid contributions to the quark self-energy [266, 267, 268]. 

One should note that the nonanalytic behaviour of (7.26) is indeed in the energy 
variable rather than the momentum, despite the fact that the derivation has been on- 
shell where the two are related, as can be verified by explicit calculation of Re /dpi, 
which is indeed analytic on the Fermi surface [270] . One consequence of this is that 
similar logarithmic terms appear in the quark-gluon vertex only in a special kinematic 
regime, namely in 

lim ; hm A^E\pi;E,p) = ^5lln J P±- + 0(\E-ri), (7.30) 

where (E' , p') and (E, p) are the four-momenta of the quarks, but not when the order 
of the limits is interchanged. Since magnetic interactions become (almost) unscreened 
only in the regime where u> -C k, it can been argued that non- Fermi- liquid corrections 
manifest themselves only through corrections on the quark self-energy rather than the 
vertices [270]. 

7.8. NLO corrections to real parts of dispersion laws at high temperature 

At nonvanishing temperature, the real parts of the dispersion laws of fcrmionic and 
gluonic quasi-particles remain IR-safe in NLO HTL-rcsummcd perturbation theory 
(in contrast to the imaginary parts which are perturbatively accessible only at T = 
or for k = and exceptional momenta where the group velocity vanishes). However, 
such calculations are exceedingly involved, and only some partial results exist so far 
in QCD [401, 402]. 

In the following, we shall restrict our attention to the case k 2 /w 2 j >• 1 and 
consider the different branches of the dispersion laws in turn. 

7.8.1. Energetic quarks and transverse gluons At large momenta only the normal 
branch of quark excitations and the transverse gluons have nonncgligible spectral 
weight. Their leading-order thermal masses become momentum independent in this 
limit and are given in one- loop order by their HTL values, m 2 ^ = rh 2 D /2 and 
Af£, = 2M 2 with rh 2 D and M 2 given by (5.18) and (5.28) respectively. While the 
HTL self-energies are no longer accurate at large momenta, their light-cone limit is 
exact to one-loop order and thus still determines the leading-order asymptotic masses 
[317, 318]. 

There are however nevertheless higher-order corrections which in fact require HTL 
resummation as indicated in (6.17) and (6.18) and which have complicated momentum 
dependence that has not been evaluated so far. Presently, only the particular averages 
(6.19) and (6.20) are known from their relation to the plasmon effect in ^-derivable 
thermodynamics as discussed in Sect. 6.3.2. 

7.8.2. Longitudinal Plasmons For momenta k 2 ^> uj 2 ^ , the longitudinal plasmon 

branch approaches the light-cone, as can be seen in figure 9. From k 2 = Usik) and 
(5.16) one finds 

w|(|k|) - k 2 (l + 4e- 6k2 /(^ 2 )-2) (7.31) 
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with e 2 = g 2 {N + Nf/2) in QCD, so the light-cone is approached exponentially as |k| 
is increased. If one also calculates the residue, one finds that this goes to zero at the 
same time, and exponentially so, too. A similar behaviour occurs in the "plasmino" 
branch of the fermion propagator at momenta k 2 3> M 2 . 

Instead of QCD, we shall consider the analytically tractable case of massless scalar 
electrodynamics as a simple toy model with at least some similarities to the vastly 
more complicated QCD case in that in both theories there are bosonic self-interactions. 
There are however no HTL vertices in scalar electrodynamics, which makes it possible 
to carry out a complete momentum-dependent NLO calculation [375]. 

Comparing HTL values of and NLO corrections to LToo = — k 2 Us/k 2 , one finds 
that as k 2 — > there are collinear singularities in both: 

nB(fc)A 2 -^ln^ (7.32) 
diverges logarithmically}), whereas 

771 

6U B /k 2 -e—= (7.33) 



(with m oc eT the thermal mass of the scalar). Because (7.33) diverges stronger than 
logarithmically, one has 5Ub > IIb eventually as k 2 — > 0. Clearly, this leads to a 
breakdown of perturbation theory in the immediate neighbourhood of the light-cone 
(/c 2 /|k| 2 < (e/ln-i) 2 ), which this time is not caused by the massless magnetostatic 
modes, but rather by the massless hard modes contained in the HTL's. 

However, a self-consistent gap equation for the scalar thermal mass implies that 
also the hard scalar modes have a thermal mass ~ eT. Including this by extending the 
resummation of the scalar thermal mass to hard internal lines renders Us/k 2 finite in 
the light-cone limit, with the result 

2 f)Q4 i 
In ' +Q(e) (7.34) 



jjresum. 

lim B 



*»-»o k 2 3k 2 
e 2 T 2 



3k 2 

such that there is a solution to the dispersion law with k 2 = at k 2 /(e 2 T 2 ) = 
|ln(2.094/e) + 0(e). Because all collinear singularities have disappeared, continuity 
implies that there are also solutions for space-like momenta k 2 < 0, so the longitudinal 
plasmon branch pierces the light-cone, having group velocity v < 1 throughout, 
though, as shown in figure 14. While at HTL level, the strong Landau damping at 
k 2 < switches on discontinuously, it now does so smoothly through an extra factor 
exp[— e-\/|k|/[8(|k| — lu)]], removing the longitudinal plasmons through over-damping 
for (|k| — u>)/|k| > e 2 . This is in fact an essential singularity in the self-energy now, 
which, as has been shown in reference [403] occurs generally at k = |k| when the hard 
degrees of freedom have nonzero mass. 

So the collinear singularities that spoil HTL-resummed perturbation theory on 
the light-cone are associated with a slight but nevertheless qualitative change of 
the spectrum of longitudinal plasmons: instead of being time-like throughout and 
existing for higher momenta, albeit with exponentially small and decreasing residue 
and effective mass, they become space-like at a particular point |k| cr ;t. ~ eTln i. For 

U This is in fact the technical reason why the longitudinal branch approaches the light-cone 
exponentially when k 2 2> w 2 ,; ■ 
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Figure 14. The longitudinal plasmon branch of scalar electrodynamics including 
NLO corrections to the HTL result. The upper of the four lines gives the HTL 
result and the lines below correspond to NLO corrections with e = 0.3, 1, and 
2, respectively. The latter three lines cross the light-cone such that the phase 
velocity starts to exceed 1, but with group velocity < 1 throughout. In the space- 
like region, the plasmon modes are damped by Landau damping, which is strong 
except in the immediat e neighbourho od of the light-cone, where it is suppressed 
by a factor of exp{— ey k/[8(k — ui)]} 



|k| > |k| cr ;t. Landau damping sets in smoothly but rapidly, so that these modes soon 
become over-damped as |k| is increased. 

This phenomenon is in fact known to occur in non-ultrarelativistic (T < m e ) 
QED [404], and has been considered in the case of QCD by Silin and Ursov [405], who 
speculated that it may lead to Cherenkov-like phenomena in the quark-gluon plasma. 

In QCD, the situation is in fact much more complicated. Under the assumption 
that the collinear singularities are removed solely by the resummation of asymptotic 
gluonic and fcrmionic thermal masses in hard internal lines, the value of |k| where 
longitudinal plasmons turn space-like has been calculated in [375]. For a pure-glue 
plasma, it reads 



Such an extended resummation can in fact be related to an improved and still gauge- 
invariant version of the HTL effective action [318]. However it is likely that damping 
effects arc of equal importance here (in contrast to scalar electrodynamics), so that 
in (7.35) only the coefficient of the logarithm is complete, but not the constant under 
the log. 

8. Resummations beyond hard thermal loops 

In the previous section we have considered perturbative corrections to dynamical 
quantities at soft scales ~ gT. As soon as there is a sensitivity to "ultrasoft" scales 
~ g 2 T, (HTL) perturbation theory breaks down and typically only leading logarithms 
can be computed as we have discussed. This is in particular the case when external 
momenta are either ultrasoft or very close to the light-cone. 



k c 2 rit . = ff 2 T 2 [ln 



1.48... 



+ ...]. 



(7.35) 
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Figure 15. Generalized ladder diagrams contributing to ultrasoft amplitudes. 



In the following we shall briefly discuss a few such cases where resummations 
of the perturbative series have to include more than the HTL diagrams and where 
important progress has been achieved recently. 

8.1. Ultrasoft amplitudes 

For ultrasoft external momenta, the infrared sensitivity may in fact become so large 
that HTL diagrams are no longer of leading order. The diagrams of figure 12 have in 
fact been evaluated in [406] for fco, |k| < g 2 T with the result 



J(v, v') In ^- — h finite 
A 



(8.1) 



*W*)--*»W"»/£/ 4 , {v . k)[v ,. k) 

with A < gT and 

/(v, v') = -*<*>(v - v') + \ 7l (V 7 ,)2 n2 . (8-2) 

7T Z y/\ - (v • V') 2 

For ho, |k| ~ g 2 T this is of the same order than the HTL self-energy ft^,, ~ rh 2 D and 
even logarithmically enhanced. 

This strong infrared sensitivity is entirely due to the HTL parts of the dressed 
vertices in figure 12, with contributions involving tree-level vertices being suppressed 
by powers of g. By blowing up the HTL vertices in figure 12, one may view the 
diagram involving one 4- vertex as one hard loop with one soft propagator insertion so 
that the largeness of (8.1) indeed means that corrections within a HTL diagram cease 
to be perturbative. This is in fact special to nonabelian gauge fields and absent in 
Abelian theories (where there are no HTL vertices involving exclusively gauge fields to 
build the diagrams of figure 12); in Abelian theories there are cancellations between 
self-energy and vertex corrections [407, 375, 408, 409], which however do not carry 
over to the nonabelian case. 

Moreover, it turns out that there are infinitely many diagrams that need to be 
taken into account in order to obtain the leading order terms in ultrasoft amplitudes. 
In the case of the gluon self-energy, such diagrams are generalized ladder diagrams as 
shown in figure 15. 

In order to characterize the leading contributions to ultrasoft amplitudes, a 
diagrammatic approach is no longer practicable, but it turns out that one can 
generalize the kinetic theory treatment of HTL's to ultrasoft amplitudes. This requires 
the inclusion of a collision term which is proportional to the colour relaxation time 
T co i ~ 1/7, where 7 is given by (7.19) in the large-momentum limit (va = 1). 

In leading-logarithmic approximation the required nonabelian Boltzmann 
equation which generalizes the nonabelian Vlasov equation (6.8) has the form 
[410, 411, 412, 413, 414, 415] 

[v ■ D, W(x, v)]° = v E a (x) + 7 J dft v ,/(v, v')W(x, v'), (8.3) 



Thermal field theory 57 




Figure 16. Ladder diagrams contributing to photon production at leading order. 
In a quark-gluon plasma, the external lines correspond to on-shell photons, the 
internal lines to dressed quark and gluon propagators. The quarks are hard but 
collinear with the external photon and exchange soft gluons such as to maintain 
collinearity. 

with I(v, v') defined in (8.2). Beyond the leading-log approximation [415, 416, 417, 
394] the quantity 7/(v, v') needs to be replaced by a collision operator that involves 
also the HTL gauge-boson self-energies. The ultrasoft vertex functions obtained 
formally in analogy to the expansion (6.9) however still share many of the remarkable 
properties of the hard thermal loops, namely gauge- fixing independence and simple 
Ward identities [415]. 

At still smaller scales, lu <C |k| <C 7, it is in fact possible to simplify the above 
effective theory of ultrasoft modes and to construct a local effective theory (Bodeker's 
effective theory) [410, 411, 418], because on distances |k| _1 much larger than the colour 
coherence length 7 -1 there are no longer propagating modes with definite colour. 
Bodeker's effective theory is a stochastic theory given by the Langevin equation 

cr c E Q = (D x B) a + C°, (C a (x)( jb (x')) = 2o- c T5 l] 5 ab 6 4 (x-x'), (8.4) 

where C lS a Gaussian noise term and o~ c the colour conductivity given to leading 
order by (7.22); beyond leading-log order, this parameter has been determined in 
reference [394]. The effective theory (8.4) is also UV finite and was put to use in the 
numerical calculation [419] of the rate of baryon number violation [122] in the hot 
symmetry-restored phase of electroweak theory, which at leading order is governed by 
nonperturbative nonabelian gauge field fluctuations with spatial momenta ~ g 2 T and 
frequencies ~ g 4 T\n(l/g). 

8.2. Light-like external momenta 

As we have seen in the example of NLO corrections to the dispersion law of longitudinal 
plasmons in section 7.8.2, the HTL perturbation theory becomes insufficient not only 
for ultrasoft momenta |k| < g 2 T, but also when amplitudes involve harder momenta 
that are however nearly light-like, which gives rise to collinear singularities. This 
problem has surfaced in particular in the calculation of the production rate of real 
(non-thermalized) photons in a quark-gluon plasma from HTL-resummed perturbation 
theory [420, 421, 422] and it turned out that damping effects on hard internal lines 
have to be included as the dominant regulator of collinear singularities [423]. 

Moreover, for a complete leading-order calculation it turns out to be necessary 
to resum all ladder diagrams built from dressed propagators of the form shown in 
figure 16. This reflects the necessity to take into account the physical phenomenon 
of Landau-Pomeranchuk-Migdal suppressionf f of photon emission [425]. However, 
contrary to the conclusion arrived at in reference [426], thermal photon production 
from a QCD plasma was shown in reference [427] not to be sensitive to ultrasoft 
nonperturbative physics, and a complete leading-order evaluation of the photon 

ft For a review in the context of energy loss calculations see e.g. reference [424] 
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production rate from a hot quark-gluon plasma was finally accomplished in references 
[428, 429]. (The analogous problem in dilepton production was recently solved in 
reference [430].) 

8.3. Transport coefficients 

Another important example where an expansion in number of loops breaks down is 
in transport coefficients, which in linear response theory are given by Kubo relations 
[431] through correlation functions in the limit of vanishing frequency and at zero 
spatial momentum of the form 



where O represents components of a conserved current or of the stress tensor. 

In scalar field theory, the required resummation (at leading order) has been 
worked out in references [432, 433, 434]. The relevant diagrams turn out to be 
uncrossed ladder diagrams, where each rung brings a pair of propagators with pinching 
singularities cut off by their thermal width. 

In gauge theory, transport coefficients turn out to depend on physics at soft 
scales such as to require also hard-thermal-loop resummation. For example, the 
leading term in the shear viscosity has the parametric form rj oc T 3 /(<? 4 ln(l/<?)), 
where the logarithm is due to the dynamical screening provided by the HTL gauge 
propagator [435, 436], and an analogous leading logarithm appears in the various 
diffusion constants [437] . The initial calculations of these quantities in kinetic theory 
however turned out to be incomplete already to leading-logarithmic accuracy and have 
been completed only recently in [438]. 

In reference [439] a purely diagrammatic rederivation of the relevant integral 
equation starting from the Kubo formula has been given, which is however incomplete 
as far as Ward identities are concerned; its completion in the Abelian case has 
been discussed in [440, 441]. To organize a diagrammatic calculation of transport 
coefficients, the 2PI approximation scheme mentioned in section 3.5 is an efficient 
means. It turns out that the 3-loop ^-derivable approximation is necessary and 
sufficient to obtain the leading-order results for shear viscosity in a scalar theory 
[442, 443], and to leading- log order in QED [443]. 

Whereas at leading-log accuracy the diagrams that need to be resummed are 
dressed uncrossed ladder diagrams like those of figure 16, the set of diagrams that 
contribute in a complete leading-order calculation is still much more complicated. 
Their resummation has nevertheless been achieved recently by means of a Boltzmann 
equation with a collision term evaluated to sufficient accuracy and with the resulting 
integral equation treated by variational techniques. The results obtained within this 
effective kinetic theory [444] so far comprise shear viscosity, electrical conductivity, 
and flavour diffusion constants [445] (bulk viscosity, however, turns out to be more 
difficult and has not yet been determined even to leading- log accuracy). 

Nonperturbative (all order in g) results for transport coefficients in a gauge theory 
have been derived for the toy model of large- flavour-number QED or QCD in reference 
[68] allowing for a test of the quality of HTL approximations, which turn out to work 
remarkably well up to the point where the renormalization-scale dependence becomes 
the dominant uncertainty. 




(8.5) 
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9. Hard thermal loops and gravity 

In ultrarelativistic field theory, the HTL approximation to the polarization tensor of 
gauge bosons determines the leading-order spectrum of gauge-field quasi-particles and 
the linear response of an ultrarelativistic plasma to external perturbations. 

In the physics of the very early universe, which is filled with a hot plasma of 
various elementary particles, the gravitational polarization tensor is also of interest. 
It describes the (linear) response to metric perturbations and its infrared behaviour 
determines the evolution of large-scale cosmological perturbations, which provide the 
seeds for structure formation that arc nowadays being studied directly through the 
anisotropics of the cosmic microwave background and which are being measured with 
stunning accuracy [446, 447]. The gravitational polarization tensor is also a central 
quantity in the theory of stochastic (semiclassical) gravity [448, 449] , which aims at 
a general self-consistent description of quantum statistical fluctuations of matter in 
a curved background geometry as a stepping-stone towards a full quantum theory of 
gravity. 



9.1. HTL gravitational polarization tensor 

If T denotes all contributions to the effective action besides the classical Einstein- 
Hilbcrt action, the energy-momentum tensor is given by the one-point 1PI vertex 
function 

r - w -7^5^ (9 ' 1) 

and the gravitational polarization tensor by the two-point function 



s (y-g(x)T^„(x) 



-U^ a0 {x,y) = Sglu>{x)5ga p iy) ~ 2 Sg^) ' (9 ' 2) 

From the last equality it is clear that n Q/3Ally describes the response of the (thermal) 
matter energy-momentum tensor to perturbations in the metric. Equating II"^ 1 ' 
to the perturbation of the Einstein tensor gives self-consistent equations for metric 
perturbations and, in particular, cosmological perturbations. 

In the ultrarelativistic limit where all bare masses can be neglected and in the 
limit of temperature much larger than spatial and temporal variations, the effective 
(thermal) action is conformally invariant, that is T[g] = T[n 2 g] (the conformal anomaly 
like other renormalization issues can be neglected in the high-temperature domain). 

This conformal invariance is crucial for the application to cosmological 
perturbations for two reasons. Firstly, it allows us to have matter in thermal 
equilibrium despite a space-time dependent metric. As long as the latter is conformally 
flat, ds 2 — (t(t, x)[dr 2 — dx 2 }, the local temperature on the curved background is 
determined by the scale factor a. Secondly, the thermal correlation functions are 
simply given by the conformal transforms of their counterparts on a flat background, 
so that ordinary momentum-space techniques can be employed for their evaluation. 

The gravitational polarization tensor in the HTL approximation has been first 
calculated fully in [450], after earlier work [451, 452, 453, 454] had attempted (in vain) 
to identify the Jeans mass (a negative Debye mass squared signalling gravitational 
instability rather than screening of static sources) from the static limit of the 
momentum-space quantity Tlfivapik), ko = 0, k — > in flat space. 
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Like the HTL self-energies of QED or QCD, n M „ Q/ 3 is an inherently nonlocal 
object. Because it is a tensor of rank 4, and the local plasma rest frame singles out 
the time direction, it has a much more complicated structure. From rf' 1 ' ', = 5® and 
fc p one can build 14 tensors to form a basis for U^uapik). Its HTL limit (ko, |k| <C T), 
however, satisfies the Ward identity 

Ak^ va0 {k) = k v T aP - k a (T aaV0v + Tp a r) av ) (9.3) 

corresponding to diffcomorphism invariance as well as a further one corresponding to 
conformal invariance, the "Weyl identity" 

rT^ap(k) = ~T aP , (9.4) 

where T a p = Pq{A5^5® — rj^ u ) and Pq the ideal-gas pressure cx T 4 . These identities 
reduce the number of independent structure functions to three, which may be chosen 
as 

n^AO^nooooWM n 2 (fc) = n / (fc)/p, n 3 (fc) = fi^(k)/ P , (9.5) 

where p = T m = 3P S b- 

The HTL limit is in fact universal: the thermal matter may be composed of any 
form of ultrarelativistic matter such as also gravitons, which are equally important as 
any other thermalized matter if the graviton background has comparable temperature. 
It reads [450] 

k° k° 4- Ikl ^ 

As in ordinary hot gauge theories, the entire HTL effective action is determined 
by the Ward identities and has been constructed in [455, 456]. The 3-graviton HTL 
vertex has moreover been worked out explicitly in [457, 458] and subleading corrections 
beyond the HTL approximation have been considered in [459]. 

From (9.6) one can formally derive an HTL-dressed graviton propagator which 
contains three independent transverse-traceless tensors. As the HTL propagators in 
gauge theories, this describes collective phenomena in the form of nontrivial quasi- 
particle dispersion laws and Landau damping (from the imaginary part of A). In 
the context of cosmological perturbations, the latter corresponds to the collisionlcss 
damping studied in [460], as we shall further discuss below. 

In the limit fcp,k 2 3> GT 4 , where G is the gravitational coupling constant, one 
can ignore any background curvature terms Gp oc GT 4 , and use the momentum-space 
propagator to read off the dispersion relations for the three branches of gravitational 
quasiparticles. One of these, the spatially transverse-traceless branch, corresponds 
to the gravitons of the vacuum theory. Denoting this branch by A, (9.6) implies an 
asymptotic thermal mass for the gravitons according to [450] 

uj\ -> k 2 + m 2 Aoc = k 2 + ^x 16nGp. (9.7) 

While in a linear response theory branch A neither couples to perturbations in 
the energy density ST 00 nor the energy flux ST 0t , the two additional branches B 
and C both couple to energy flux, but only C to energy density. The additional 
branches correspond to excitations that are purely collective phenomena, and like 
the longitudinal plasmon in gauge theories, they disappear from the spectrum as 
k 2 /(GT 4 ) — ► oo. In contrast to ordinary gauge theories, however, they do so by 
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acquiring effective thermal masses that grow without bound, while at the same time 
the residues of the corresponding poles disappear in a power-law behaviour (instead 
of becoming massless and disappearing exponentially) [450]. (Some subleading 
corrections to these asymptotic masses at one-loop order have been determined in 
[461], but these are gauge dependent and therefore evidently incomplete.) 

For long wavelengths k 2 ~ GT 4 the momentum-space HTL graviton propagator 
exhibits an instability, most prominently in mode C, which is reminiscent of the 
gravitational Jeans instability and which is only to be expected since gravitational 
sources unlike charges cannot be screened. However, this instability occurs at 
momentum scales which are comparable with the necessarily nonvanishing curvature 
R ~ GT 4 . The flat-space graviton propagator is therefore no longer the appropriate 
quantity to analyse; metric perturbations have to be studied in a curved time- 
dependent background. 



9. 2. Self consistent cosmological perturbations from thermal field theory 

The conformal covariance of the HTL gravitational polarization tensor IT MJ , Q/ 3(x, y) as 
expressed by the Weyl identity (9.4) determines this non-local quantity in a curved 
but conformally flat space according to 

n M1/ce/ 3(x,2/) = o-(x)o-(y)n^ al 3(x - y) . (9.8) 

A closed set of equations for metric perturbations is obtained from using this 
in the right-hand side of the Einstein equation linearized around a conformally-flat 
background cosmological model [450, 462, 463] 

$9 p Jx< 5g a p{x') 

= AixG [T^g a ? + 2T" V] Sg a p(x) + 4== / W va f*{x, x')5g af} (x>). (9.9) 

V-9{x) Jx< 

In a hydrodynamic approach, 8T^ V is usually determined by certain equations of state 
together with covariant conservation, the simplest case of which is that of a perfect 
fluid, which has been studied in the pioneering work of Lifshitz [464, 465]. Many 
generalizations have since been worked out and cast into a gauge invariant form by 
Bardeen [466] (see also [467, 468]). 

Relativistic and (nearly) collisionless matter, however, has more complicated 
gravitational interactions than a perfect fluid. This is usually studied using classical 
kinetic theory [469, 470] and the case of purely collisonless matter has been worked 
out in [471] with some numerical solutions obtained in [472, 460, 473] for particular 
gauge choices. 



9.2.1. Purely collisionless matter In [462], the self-consistent equations for (scalar) 
density perturbations of a radiation-dominated Robertson- Walker-Friedmann model 
with collisionless matter have been derived using the HTL gravitational polarization 
tensor. By virtue of the diffeomorphism Ward identity (9.3), these equations turn out 
to be automatically gauge independent. 

For example, in the simple radiation-dominated and spatially-flat Einstein-de 
Sitter model 

ds 2 = a(T)(dr 2 - dx 2 ), <t(t) = ( 9 .10) 
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where p is the energy density when a = 1 and r is the conformal time (which equals 
the size of the Hubble horizon in comoving coordinates), the scalar part of metric 
perturbations can be parametrized in terms of four scalar functions 



' V. Dj ASij + I',.,, , 

Of these, two can be gauged away by diffeomorphisms. But in a gauge- 
invariant framework only gauge-invariant combinations enter nontrivially. Only two 
independent gauge-invariant combinations exist, which may be chosen as 

$ = A + ^( D - \B) (9.12) 

II = \{B + —B + C — A) — D — —D, (9.13) 
a a 

where a dot denotes differentiation with respect to the conformal time variable r. 

Each spatial Fourier mode with wave vector k is related to perturbations in the 

energy density and anisotropic pressure according to 

5 = ^x 2 $, 7r anis . = ix 2 II, (9.14) 

where 

* skr =>7<ty (915 > 

which measures the (growing) size of the Hubble horizon over the wavelength of a 
given mode (which is constant in comoving coordinates). In (9.14) energy density 
perturbations S are defined with respect to space-like hypersurfaces representing 
everywhere the local rest frame of the full energy-momentum tensor, whereas 7r an i s . is 
an unambiguous quantity, since there is no anisotropic pressure in the background. 

Correspondingly, when specifying to scalar perturbations, there are just two 
independent equations contained in (9.9). Because of conformal invariance, the trace 
of (9.9) is particularly simple and yields a finite-order differential equation in x, 

$" + -$' + -$ = -n - -IT (9.16) 

a: 3 3 a; 

(a prime denotes differentiation with respect to the dimensionless time variable x). The 
other components, however, involve the nonlocalities of the gravitational polarization 
tensor. These lead to an integro-differential equation, which upon imposing retarded 
boundary conditions reads [462] 

{x 2 - 3)$ + 3x& = 611 - 12 / dx j (x - x')[&{x') + IL'(x')] + ip{x - x ) (9.17) 

where jo(x) = sin(x)/x arises as Fourier transform of Ui(oj/k) in (9.6). <p(x — x ) 
encodes the initial conditions, the simplest choice of which corresponds to (p(x — xq) cx 

j (x - Xq). 

Similar integro-differential equations have been obtained from coupled Einstcin- 
Vlasov equations in particular gauges, and the above one can be shown to arise from a 
gauge-invariant reformulation of classical kinetic theory [474] . Usually, such equations 
are studied numerically, with only some asymptotic behaviour having been analysed 
analytically [475, 476]. Remarkably enough, they can be solved analytically [462] 
provided initial conditions are formulated for x — > 0. In this case a power series 
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Figure 17. The energy-density contrast (arbitrary normalization) as a function 
of x/tt for a collisionless ultrarclativistic plasma (full line), a scalar plasma with 
quartic self-interactions \<fi 4 and A = 1 (dashed line), and a perfect radiation fluid 
(dotted line). 

ansatz for $ and II leads to solvable recursion relations for an alternating series that 
converges faster than trigonometric functions. 

This also holds true for the vector (rotational) and tensor perturbations and when 
the more realistic case of a two-component system of a perfect radiation fluid combined 
with a collisionless ultrarelativistic plasma is considered [463] . In the case of rotational 
perturbations these studies led to novel solutions not considered before [477] . 

In figure 17, the solution for the energy-density contrast is given in a doubly- 
logarithmic plot (full line) and compared with the perfect-fluid case (dotted line). In 
the latter, one has growth of the energy-density contrast as long as the wavelength 
of the perturbation exceeds the size of the Hubble radius (x <C 1)- After the Hubble 
horizon has grown such as to encompass about one half wavelength (x — n) , further 
growth of the perturbation is stopped by the strong radiation pressure, turning it 
into an (undamped) acoustic wave propagating with the speed of sound in radiation, 
v = 1/V3- The collisionless case is similar as concerns the superhorizon-sized 
perturbations, but after horizon crossing, there is strong damping ~ 1/x, and the 
phase velocity is about 1. This indeed reproduces the findings of the numerical studies 
of reference [460]. They can be understood as follows: a energy-density perturbation 
consisting of collisionless particles propagates with the speed of their constituents, 
which in the ultrarelativistic case is the speed of light, and there is collisionless 
damping in the form of directional dispersion. 

Tensor perturbations, which correspond to primordial gravitational waves, are 
also modified by a nearly collisionless background component as worked out in 
references [463, 474]. 
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(a) 

Figure 18. Two examples of infrared divergent contributions to the gravitational 
polarization tensor in scalar <f> 4 theory beyond two-loop order. 

9.2.2. Weak self-interactions and HTL resummation The thermal-field-theory treat- 
ment of the effects of an ultrarelativistic plasma on the evolution of cosmological 
perturbations can be used also to study the effect of weak self-interactions in the 
plasma by calculating higher-order contributions to the gravitational polarization 
tensor. A virtue of this approach is that everything is formulated in purely geometrical 
terms, without explicit recourse to perturbations in the (gauge variant) distribution 
functions of a kinetic-theory treatment. 

In reference [63], the gravitational polarization tensor has been calculated in a 
\<j) 4 theory through order A 3 / 2 . The next-to-leading order contributions to Tl^af) at 
order A 1 are contained in the high-temperature limit of two-loop diagrams and their 
evaluation is straightforward. However, starting at three-loop order, there are infrared 
divergences which signal a breakdown of the conventional perturbative series. This 
is caused by the generation of a thermal mass oc \/~\T for the hot scalars. If this is 
not resummed into a correspondingly massive scalar propagator, repeated insertions 
of scalar self-energy diagrams in a scalar line produces arbitrarily high powers of 
massless scalar propagators all with the same momentum, and thus increasingly 
singular infrared behaviour (figure 18a). 

However, it is not sufficient to resum this thermal mass for the hot scalars, as this 
would break conformal invariance (which in the ultrarelativistic case is only broken 
by the nonthermal conformal anomaly). Indeed, there arc also vertex subdiagrams 
oc AT 2 that have a similar effect as a self-energy insertion, see figure 18b. As in 
the hard-thermal-loop resummation program for ordinary gauge theories, one has to 
resum also nonlocal vertex contributions. Doing so, the result turns out to satisfy 
both the diffeomorphism and conformal Ward identities. 

In the low-momentum limit of interest in the theory of cosmological perturbations, 
the function Hi in (9.6) that governs the evolution of scalar perturbations reads, to 
order A 3 / 2 

1 5 5A I" / 1x2 1 " 2 

Hi(zu) = ?u artanh + — — ■ 

zu 4 air 2 



. i v ,1 ^ 

2 I zu artanh — — zu artanh — 

zu I zu zu z - 



5^3/2 r / | 

+ - 3 3 [zu 2 — 1 — zu V zu 2 — lj izu artanh — 



+6 \zu\J zu 2 — 1 — zu 2 | zu artanh — 

V Vcc 2 - 1/ 07 

2 

ZU ZU ZU / 2 

+ T~2 TW2 + 3— 2 7 + 6 / , = - SzuVzu 2 - 1 + 3tu 

{ZU 2 - I) 6 / 2 ZU 2 - I V^7 2 - 1 



(9.18) 



where zu = fc /|k|, and similarly complicated expressions arise for n 2 and n 3 , which 
in the collisionlcss limit (9.6) are simple constants. 

The Fourier transform of this expression determines the kernel in the convolution 
integral of (9.17). At order A 1 , it can still be expressed in terms of well-known special 
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functions [478], whereas at order A 3 / 2 this would involve rather intractable integrals 
over Lommel functions. However, all that is needed for finding analytical solutions 
is their power series representations which are comparatively simple. Given them, it 
is as easy as before to solve the perturbation equations. However, one finds that the 
asymptotic behaviour x 3> 1 eventually becomes sensitive to higher and higher loop 
orders. The reason for this is that higher loop orders come with increasingly singular 
contributions at w = ±1 to H(ru), and the large- a; behaviour is dominated by the 
latter. This could be cured by a further resummation similar to the one considered 
in the case of the dispersion relations of longitudinal plasmons in the vicinity of the 
light-cone in (7.34), but it turns out that a particular Pade-approximant based on the 
pcrturbative result reflects the effects of this further resummation quite well [479] . The 
result for the density perturbations in a scalar plasma with A</> 4 -interactions and A = 1 
are shown in figure 17 by the dashed line, where it is compared with the collisionlcss 
case (full line) and the one of a perfect radiation fluid (dotted line). The effects of the 
self-interactions within the ultrarelativistic plasma become important only for x > n, 
where the strong collisionless damping is somewhat reduced and the phase velocity is 
smaller than 1. 

10. Conclusion 

In this review we have concentrated on the progress made during the last decade 
in calculating by analytic means static equilibrium and dynamic near-equilibrium 
properties in ultrarelativistic gauge theories. Even when the practical utility of 
the obtained results is often open to debate, they are hopefully paving the way for 
understanding the more complicated realistic situations to be analysed in present or 
future heavy-ion experiments, or in astrophysical problems. 

Static quantities at high temperature and not too high chemical potential can be 
efficiently investigated by lattice simulations, whereas perturbative methods seemed to 
be too poorly convergent to be of any predictive power in the applications of interest. 
However, in the last few years there have been several different investigations leading 
to the conclusion that a careful resummation which emphasizes a weakly interacting 
quasiparticle picture, or, in dimensional reduction, the cffectivc-ficld-theory aspect, 
is able to provide dramatic improvements, restoring prcdictivity even down to a few 
times the deconfinement temperature in strongly interacting QCD. 

In nonabelian gauge theories the perturbative approach is limited by the magnetic 
mass scale, or, more precisely, the physics of confinement in the chromo-magnetostatic 
sector. However, weak-coupling effective-field-theory methods can still be used to 
combine analytical and (numerical) nonperturbative techniques to achieve further 
progress which is complementary to or perhaps beyond the capacities of a direct 
numerical approach in four dimensions. 

In the case of dynamic properties, and also in the case of static properties 
when dimensional reduction is not applicable as in cold ultradegenerate plasmas, the 
analytical approach is even more important, and it is in fact here that the greatest 
variety of phenomena are encountered. At weak coupling, there are several spatial 
or temporal scales that need to be distinguished and are at the root of the required 
resummations in a perturbative treatment. In terms of frequencies and momenta, the 
first important scale below the hard scale of temperature or chemical potential is the 
soft scale set by the Debye mass responsible for screening of electric fields as well as 
for the frequency of long- wavelength plasma oscillations. This is the realm of HTL 
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resummations, which is however limited by an eventual infrared sensitivity to ultrasoft 
scales occurring at some (mostly very low) order of the expansion. 

Resummations at zero temperature and high chemical potential are not limited by 
the magnetic mass scale, but nearly static magnetic modes lead to qualitative changes 
(non- Fermi-liquid behaviour) . There is furthermore the nonpcrturbative phenomenon 
of colour superconductivity, which is however to some extent accessible by weak- 
coupling methods, which strongly depend on HDL resummation techniques. 

At high temperature, there is, between soft and ultrasoft scales, a further scale set 
by the damping rate of the hard plasma constituents, which in a nonabelian plasma 
determines a colour coherence length. The corresponding energy scale is enhanced 
over the ultrasoft scale by a logarithm in the inverse coupling, which allows for novel 
systematic developments that we could only cursorily describe. 

In the final brief excursion to general relativity, we described the role of the HTL 
contributions to the gravitational polarization tensor in the theory of cosmological 
perturbations. There the soft scale is given by the Jeans mass which is in a 
cosmological situation comparable with the scale of the (inverse) Hubble horizon. 

It should be needless to emphasize that many interesting topics in thermal field 
theory have been covered only cursorily or even not at all. Some of the very recent 
developments in fact are about to leave the arena of traditional thermal field theory 
towards a more general, fully noncquilibrium field theory, which is only timely in 
view of the wealth of experimental of data to be expected from the modern relativistic 
heavy-ion colliders. But (near-) equilibrium thermal field theory will certainly continue 
to be an important theoretical laboratory where many physical concepts are brought 
together for cross-fertilization. 



A. Spectral representation of HTL/HDL propagators 

A. 1. Gauge boson propagator 

For the two nontrivial structure functions of the HTL/HDL gauge boson propagator 
corresponding to the branches A and B it is convenient to separate off the kinematical 
pole at k 2 = in A# and to define 

A = -A - = r dk '° P t{k ' M (A 1) 

A fcg - k 2 - II A (fco, |k|) J-oo 27T k' -k ' 

A =_^!a -1 _ [°° dk' p e (K,\k\) 1 

" k2 + H B (fc ,|k|) loo 27T k' -k k 2 ' {A - Z) 

where IL4,b are the HTL quantities given in (5.15), (5.16). The spectral functions are 
given by 

p t .e(k , |k|) = DiscA M (fc ,|k|) = 2 lim Im A M (fc + ie, |k|) 

= 27re(fc ) M (|k|) 5(k 2 - ^(|k|)) + /3 M (fc , |k|) 6(-k 2 ) (A.3) 

with w 2 ^(|k|) as shown in figure 9 (for k 2 > 0). For small and large values of k 2 they 
are approximated by [217] 

^t-^li. + j|k 2 , w?^w= L + ^k 2 , k 2 «w 2 1 . (A.4) 

u? t ~ k 2 + ml,,, ujj ~ k 2 + 4k 2 exp (~\- - 2 V k 2 > (A.5) 
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where uu 2 x = m 2 D /3 is the plasma frequency common to both modes, and ra 2 ^ = rh 2 D /2 
is the asymptotic mass of transverse quasiparticles. The effective thermal mass of 
mode I (or B) vanishes exponentially for large k 2 . 
The residues z t j, are defined by 



z t.t 



d_ 



and explicitly read 

z t = 



k a =u t (|k|) 



(A.6) 
(A.7) 



m 2 D k 2 - (fc 2 ) 2 fe = Wt (|k|)' k 2 (rh 2 D -k 2 ) 
with the following asymptotic limits [218] 
4k 2 w 2 , 

/, 4k 2 „\ 4k 2 , * . \ , , 




^-^r^- 2 ;- ^^ exp lr^- 2 >' k>>m - (A - 9) 

The singular behaviour of for k 2 — > is in fact only due to the factor fc 2 /k 2 in 
(A. 2); the residue in Ab approaches 1 in this limit. For k 2 ^> m 2 ^, the residue in 
vanishes exponentially, as mentioned in Section 5.1.2. 
The Landau-damping functions (3 t ,e are given by 

(3 t (k , |k|) = nm 2 D ^p-\A t (k , |k|)| 2 , fa(k , |k|) = nm 2 D ^\A e (k , |k|)| 2 . (A.10) 

These are odd functions in fco which vanish at fco = and at k 2 , = k 2 . For large k 2 
and fixed ratio fco/|k|, f3 t .i decay like 1/k 4 . 

The spectral functions pt.t satisfy certain sum rules which can be obtained by 
a Taylor expansion of (A.l) and (A. 2) in fco [307, 480, 10]. A special, particularly 
important case is obtained by putting fc = in (A.l) and (A. 2) which yields 



(A.11) 



r dk Q p t (k ,\k\) _ 1 r° dfco p/(fc ,|k|) _ 1 1_ 

./-oc 2tt fc k 2 ' 2ir fc k 2 k 2 +rh 2 D 

More complicated sum rules have been found in applications of HTL 
resummations. In [481] it has been shown that one can reexpress the integrals 
involving only the continuous parts of the spectral functions appearing in the energy 
loss formulae of heavy particles [482] in terms of generalized (Lorentz-transformed) 
Kramers-Kronig relations; in [483] a sum rule involving also only the Landau damping 
domain \x\ < 1 with x = fco/|k| has been derived which appears in calculations of the 
photon or dilepton production rates in a quark-gluon plasma [484] 

1 f 1 dx 2imrL(x) 1 1 



(A.12) 

[t + ReIL»] 2 + [ImIL(x)] 2 t + RcIL(oo) t + ReIL(0) 
for t > and i = A, B. 

A peculiar sum rule has been encountered in [99], 

2Imn A ReA t - Imn B ReA4 = (A.13) 

which has so far only been shown to hold by numerical integrations. Like (A.12) this 
sum rule receives contributions only from the Landau damping domain fc 2 < 0, but it 
involves the two branches at the same time and it holds only under both fc and |k| 
integrations. 
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A. 2. Fermion propagator 

The spectral representation of the two branches of the HTL fermion propagator is 
given by 



fc T(|k|+S±(fc ,|k|) J-oo^ k' -k 
where p± are defined in analogy to (A. 3), 

p±{ko, |k|) - 27re(fc ) z ± (\k\) S(k 2 - w|(|k|)) + P±(k , |k|) fl(-fc 2 ) (A.15) 

with w±(|k|) as shown in figure 10. For small and large values of k 2 they are 
approximated by [228] 

w+~M+^, w_~M-^, k 2 <M 2 (A.16) 

c 2 ~k 2 + M^, w _~|k|fl + 2expf-^--l)y k 2 »M^ (A.17) 



M 2 

where = 2M 2 is the asymptotic mass of energetic fermions of the (+)-branch. 
The effective thermal mass of the additional (— )-branch vanishes exponentially for 
large k 2 . 

The residues z± are given by the simple expression 

z ± = (A.18) 
2M 2 V ^ 

with the following asymptotic limits [218] 

^■4( k l4 --"§ exp (-€-')• k2>>AC <A ' 20) 

For k 2 3> M^, the residue in A_ vanishes exponentially. 
The Landau-damping functions (3± are given by 

/3± =ttM 2 M^| A± (fc ,|k|)| 2 . (A.21) 

For large k 2 and fixed ratio fco/|k|, they decay like 1/k 2 . 

In contrast to the gauge boson case, the spectral functions are not odd functions 
in ko but rather obey 

p + (-fc ,|k|)=p_(fc ,|k|). (A.22) 

Sum rules for these spectral functions have been discussed in detail in [10]. A more 
complicated one that plays a role in the HTL resummed calculation of the hard photon 
production rate [485, 486] has been given recently in [481]. 
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